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A STUDY OF THE MICROSTRUCTURE AND 
OPTICAL PROPERTIES OF THIN LEAD-DIELECTRIC 
CERMET FILMS 


I. INTRODUCTION 


Thin cermet films, similar to other vacuum-deposited thin films, are formed 
through a process of nucleation that starts with the appearance of isolated 
three-dimensional nuclei [1-5]. If deposition continues, the nuclei grow, intergrow, and 
finally form a continuous deposit film at an average thickness of 5 to 100 nm [6] . 

Previous studies [7-9] have shown that composite dielectric-metal (cermet) films 
consist of metallic grains embedded in a dielectric network. Miller and Shim [7] studied 
the conductive properties of co-sputtered Au-Si0 2 , thin film cermets. The electric 
resistance was found to vary exponentially with gold concentration. The microstructure 
was examined briefly by transmission electron microscopy and was found to consist of 
discontinuous networks of Au and Si0 2 . A similar study of flash-evaporated Cr-SiO films 
involving film structure, resistivity, and temperature coefficient of resistance was carried 
out by Gland, et al. [8] . Electron microscopy revealed a discontinuous microstmcture that 
varied with metallic content. Another study of the structure and resistivity of co-sputtered 
metal-dielectric films was made by Beckerman and Thun [9]. Various metal-dielectric 
combinations of percentages and materials revealed a discontinuous, grainy microstructure. 

Because of their microstructure, cermet films are intimately related to the 
extensively studied but poorly understood ultrathin metallic films which consist entirely of 
the initial, isolated nuclei or grains. Therefore, we shall briefly review the studies made of 
ultrathin metal films, most of which are concerned with the electrical conductivity of the 
films. 


The problem of conduction by tunneling between isolated, conducting grains was 
analyzed in 1930 by Frenkel [10], Further studies on conduction through thin dielectric 
barriers were made by Emtage and Tantrapom [11], who considered Schottky emission 
through thin insulating films; Stratton [12], who extended previous tunneling theory by 
including temperature dependence and image force corrections; Hartman, et al. [13] , who 
made extensive experimental measurements of electrical conduction through SiO films; Hill 
[14], who did an excellent survey of the various proposed mechanisms and relevant 
experimental results; and Gundlach and Simmons [15] , who analyzed the range of validity 
of the WKB approximations used in some of the proposed mechanisms. 

The first detailed investigation of conduction in ultrathin metal films was done by 
Neugebauer and Webb [16] , who concluded that the conduction process consisted of, first, 
thermally activated charge carrier creation involving charge transfer between initially 
neutral particles and, second, the drift velocity of those charges in an applied field. Charge 


transfer between charges occurred by tunneling. Neugebauer’s and Webb’s mechanism was 
supported by experimental results. An alternate mechanism was then proposed by Hartman 
[17], who assumed tunneling between neutral grains. Hartman’s mechanism, which was 
also in agreement with experiment, involved an activation energy corresponding to 
differences in energy levels of the band structure of the individual particles, which was 
assumed to be discrete because of the small size of the particles. The Neugebauer and Webb 
mechanism suggested that a certain number of the microparticles must be charged, thus 
requiring a thermal activation energy to move the appropriate electrons, with tunneling 
then occurring between charged and uncharged microparticles. Workers agreed that some 
kind of activation energy was required to reconcile theory and experiment, but it was 
questionable as to which theoretical form it would take. Experimental results slightly 
favoring the Hartman theory were obtained by measurements of electrical resistance-strain 
characteristics of thin metal ‘films made by Parker and Krinsky [18], but then a further 
modification of the Neugebauer and Webb thermal activation mechanism by Weitzenkamp 
and Bashara [19] brought that theory more in line with experiment. 

At this time a different mechanism was proposed by Herman and Rhodin [20], 
who suggested that conduction should include contributions caused by electron tunneling 
between regions in the substrate surface immediately under the metallic particles. Charge 
carriers would be contributed to the insulator surface regions by the metallic 
microparticles. Their activation energy, which also produced results agreeable with 
experiment, was electrostatic in nature and depended upon the average size and separation 
of the metal grains. 

A modification of the tunneling process common to all the previously described 
mechanisms was made by Milgram and Lu [21 ] , who assumed that trapping centers existed 
in the substrate. They considered charge carriers created both by thermal excitation and by 
field injection coupled with a charge transportation mechanism consisting of tunneling via 
the trapping centers. 

Further modification of the Neugebauer and Webb model was done by Swanson, et 
al. [22], who considered electron transfer between charged particles other than 
charged-to-neutral particle transfers. 

The confusion was finally resolved by Hill [23], who, in an extensive theoretical 
and experimental work, established a mechanism similar to that of Neugebauer and Webb, 
but with an activation energy that considered substrate electrical properties and 
microparticle dimension and a tunneling mechanism that incorporated temperature 
dependence and the effects of an active substrate. 

It should be noted that an interesting proposal by Morris [24] suggested correlated 
conductivity and Hall-effect measurements on cermet films as a method of finally 
determining which of the above mechanisms was correct. Unfortunately cermet films never 
showed a measureable Hall effect. 1 


1. J.E. Morris, Private Communication, April 28, 1970. 
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The only examinations of cermet films have been in the form of conduction and 
structural studies. In addition to those studies mentioned, a detailed conduction study was 
made by Christen and Hewitt [25] , who measured resistivities and temperature coefficients 
of resistance for radiofrequency, co-sputtered gold-silicon dioxide and vacuum-evaporated 
gold-silicon monoxide films. Their results are compatible with the mechanism of 
Neugebauer and Webb. 

No optical studies have been made on thin cermet films until the present study, 
although such films would be roughly included in the general study of the optical 
properties of thin absorbing films made by Abeles [26] . Related optical studies have been 
made by workers examining ultrathin metal films. 

The theory of Maxwell Garnett [27] applies to island structure films, such as 
ultrathin metal films. Studies of the optical properties of such films have been made and 
are reviewed by Heavens [28] . Maxwell Garnett assumed spherical metal particles that were 
small with respect to a wavelength of light. His results agree with the more general theory 
of Mie [29] for light absorption and scattering by an arbitrarily -sized sphere in the limit of 
small, widely-separated spheres. Some investigators [30, 31] reported disagreement with 
the Maxwell Garnett theory and assumed nonspherical microparticles without any direct 
evidence for the necessity of this assumption. The resultant treatments are extremely 
complex. However, it was shown in the later work of Doremus [32, 33] that previous 
results could be understood in terms of the Maxwell Garnett theory and that it was not 
necessary to assume that the particles were nonspherical or to make any other 
modifications of the theory. It might be noted that while in the original formulation of the 
theory it was stated that for thin films the Maxwell Garnett expression for the dielectric 
constant might have to be modified for light falling on the film at an oblique angle, the 
results of Doremus indicate that no such modification is necessary. It is shown in the 
present study that the Maxwell Garnett model applies with equal validity to thin cermet 
films. 


In this study, variations in surface structure are correlated with variations in the 
internal metallic distribution. This is done by electron microscopy, using both replicating 
techniques and direct observation of films, and by microprobe analysis. The structural 
model found is approximated by the Maxwell Garnett model, which is then used to 
generate theoretical refractive indices. Rough indices of refraction experimentally using the 
Tekucheva method [34] for interpretation of transmission and reflectivity measurements 
and more exact experimental values found using ellipsometry [35, 36] are compared to the 
theoretical Maxwell Garnett refractive index values. 


II. MICROSTRUCTURE EXAMINATION 


This section presents the experimental procedure used in this study for the 
examination of the microstructure of thin cermet films. 
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A. Film Formation 


The films were made in a Consolidated Vacuum Corporation (CVC) Plasmavac 
AST-200 sputtering unit, which may be seen in Figure 1 . The unit utilizes an RF-sputtering 
technique which enables one to sputter dielectrics as well as conductors. This technique is 
well discussed in Holland [37]. Plasmavac chamber pressures were less than 10~ s torr 
before argon was introduced and sputtering began. 



fc 

4 » 

0 0 ] 




Figure 1 . Photograph of CVC Plasmavac AST-200 
sputtering unit. 

polished before each deposition to remove any lead oxide. 


Single sputtering tar- 
gets were used. This tech- 
nique should ultimately yield 
films more uniform and 
reproducible than those 
resulting from the dual-target 
techniques used by others [7, 
8] . The targets consisted of a 
nickel backing, lead glass, and 
varying amounts of pure lead. 
A typical target was formed 
by firing a coating of Corning 
7572 or 7575 lead glass onto 
a 0.157-cm-thick square of 
nickel and fusing strips or 
dots of pure lead to the glass. 
The glass fired at about 
500° C. This technique was 
suggested by R.F. DeHaye, 
Marshall Space Flight Center 
(MSFC). Three targets were 
made, having 10-, 25-, and 
50-percent pure lead by sur- 
face area, respectively, 
soldered to their fronts. The 
targets were mechanically 


Film thickness and sputtering rates were monitored by a Sloan 103-850 quartz 
crystal monitor. This unit uses a newly developed round quartz crystal which mounts inside 
a shielded stainless steel cap. Both the crystal mounting surface and its signal contact are 
shielded from the evaporant stream. Thus, unlike previous monitors, the unit gives stable 
readings while in a plasma. Monitor calibration was done by Tolansky interferometry [38] 
in a Sloan M-100 Angstrometer. Measurement by quartz crystal monitor is described by 
Holland [37] . Sputtering rates of about 10 nm/min were typical. 

Films with thicknesses of 20, 40, 60, and 80 nm were made from each of the 
targets, and additional films of 100 nm and 100 and 140 nm were made from 25- and 
10-percent lead targets, respectively. 
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B. Examination by Microprobe 


An Advanced Metals Research (AMR) Corporation X-ray spectrometer attachment 
for the Philips EM-300 transmission electron microscope was used to measure the actual 
lead content of the films. All films measured were deposited onto electron microscope grids 
that were precovered with a Formvar film. Thus, the films could be inserted directly into 
the electron microscope for analysis. A 40 nm-thick, evaporated-lead film was used as a 
standard, and a simple ratio measurement of peak heights was used to determine the 
percentage of lead in the cermet films. Peak heights are taken to be proportional to lead 
content. This procedure is described in References 39 and 40. 

The AMR X-ray spectrometer attach- 
ment and the associated Norelco Mark III 
Data Control and Processor may be seen in 
Figure 2. With this equipment attached to 
the electron microscope, one bombards the 
sample with the microscope’s electron beam 
and examines the intensity of the emitted X- 
rays as a function of the X-ray energy. 
Because of the discrete emissions character- 
istic of any given element, one is able to 
determine which elements are present in the 
sample by comparing the resultant peaks or 
lines with known, tabulated values. By com- 
paring the intensities of the lines with line 
intensities from the appropriate known 
standards, it is also possible to determine 
quantitative values for sample composition. 


C. Examination by Electron Microscopy 

Replicas suitable for electron micro- 
Figure 2. Photograph of AMR X-ray scopy were made from each film. This was 
Spectrometer attachment and Norelco done by the following well-established proce- 
Mark III Data Control and Processor. dure [41], Presoftened replicating tape was 

pressed onto each film and left to set for at 
least 3 hours. These surface impressions were then removed, coated with a thin 
vacuum-sputtered carbon film which followed the contours of the impressions, and 
shadowed with chromium from a source angled at 30 deg with respect to the replica plane. 
The replicating tape was then desolved away in an acetone bath, leaving the shadowed 
carbon replicas which were scooped out of the bath onto a microscope grid. The resulting 
specimens show the surface structure of the films under study. 

It was also possible to examine the films directly by electron microscopy. Some 
films were deposited directly onto 3.05-mm, 200-mesh Fullum electron microscope grids 
that were precovered with a Formvar film. These specimens could be inserted directly into 
the electron microscope and examined. All electron microscope examinations were made 
on a Philips EM-300 transmission scope, which may be seen in Figure 3. 
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Figure 3. Photograph of Phillips EM-300 transmission electron microscope. 

III. MEASUREMENT OF OPTICAL PROPERTIES 

This section presents the experimental techniques used in this study for 
measurement of optical properties of thin cermet films. The films were formed in the same 
manner as before except that they were deposited onto Kodak B351 cover glass substrates. 

A. Through Absorption and Reflection 

It is possible to determine thin Film indices of refraction by measuring the intensity 
of monochromatic light reflected by and transmitted through both the bare substrate and 
the film-substrate system. The approach used in this study is the method of Tekucheva 
[34] , whose method requires knowledge of the light source wavelength, the film thickness, 
the substrate index of refraction, and the appropriate transmission and reflection 
measurements. For this study, Kodak B351 slide cover glass was used as substrates. The 
index of refraction of Kodak B351 glass has been determined to be 1.538 ± 0.001 at 546.1 
nm, 2 and this value as well as the film thicknesses measured by the Sloan 103-850 quartz 
crystal monitor were the values utilized in the calculations in this study. These were also 
the values used in evaluating ellipsometric measurements. 

2. J.L. Zurasky, Private Communication, Marshall Space Flight Center, Alabama, 1970. 
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A Cary Model 14 recording spectrophotometer was used to make transmission and 
reflection measurements. The instrument automatically records transmission and reflection 
spectra over the 186 to 2600 nm wavelength region with accuracies of ±3 percent,- using 
assorted light sources, energy receivers, slits, collimators, prisms, and gratings. In this 
report, measurements were confined to the 325- to 800-nm region, where a tungsten lamp 
acted as source and a 1P28 multiplier phototube acted as energy receiver. 

Measurements were taken on all the cermet films under study, but meaningful 
results were obtained only for the more absorbent films from the 50-percent pure lead 
sputtering target, since the accuracy of the Cary Model 14 was insufficient to detect the 
transmission variances required to use the Tekucheva method. A presentation of this 
method follows. However, the values obtained acted as a rough guide in the determination 
of the iteration range utilized in the McCrackin program for solution of the Drude 
equations. 

Refractive Indices by the Tekucheva Method. Tekucheva ’s [34] method allows one 
to obtain refractive indices of thin films from measurements of film transmission and 
reflection. A derivation of this method follows. 

Consider an absorbing medium of thickness d, reacting with monochromatic light 
of wavelength X 0 and initial intensity I 0 . We know that the transmitted light has 
intensity I, where 


I 


Io 


-ad 

e 


( 1 ) 


and a is the absorption coefficient of the medium. We may relate a to the medium’s 
optical constants by 


a = 47r n/c/Xo 


( 2 ) 


where 


n = refractive index of medium 
k = absorption index of medium. 

Note that the complex refractive index n’ would be defined here as 


n* = n ( 1 - ix) 


(3) 
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Equation (1) is valid only if there is little reflection at the boundaries. For a medium 
without a substrate but with reflection at both boundaries, equation ( 1 ) becomes 

I - I 0 ( 1 - r) 2 e -Q?d , (4) 


where r is the reflection coef- 
ficient for the air-medium 
boundary, and the incident 
light is normal. 

Now consider an 
absorbing film on a substrate, 
as shown in Figure 4. Here, we 
have reflection at three bound- 
aries — the air-film, the film- 
substrate, and the substrate-air 
boundaries. Also, in addition 
to film absorption we have 
absorption by the substrate. 

Assuming film and substrate 
thicknesses of d! and d 2 
and film and substrate absorp- 
tion coefficients of a! and 
a 2 j respectively, we may 
write 

Ii = I 0 U - RiMl - R 2 m - R 3 )e"" ldl e _0;2d2 , (5) 


where Ri , R 2 , and R 3 are the reflection coefficients at boundaries 1, 2, and 3, 
respectively, as shown in Figure 4, and Ii is the intensity of the light transmitted by the 
film and substrate. We are considering the incident light to be normal, and we are 
neglecting multiple-reflection terms, which are relatively small at normal incidence for the 
films under study. 

We note that if we consider the intensity of the light transmitted by the bare 
substrate only, I 2 , we may write 

I 2 = I 0 (l - R 3 ) 2 e -Q!2d2 . (6) 



Figure 4. Film substrate system in unprimed mode. 
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Solving for e _Q!2 ^ 2 in equation (6) and substituting into equation (5) yield 


It - I 2 (l - Ri)(l - R 2 )e"“ lCll /(l - R 3 ) , 


-«i d i Ii( l - r 8 ) 


I 2 (l - R t )(i - r 2 ) ' 


We define the overall reflection coefficient for reflection from all three boundaries as 


R = i r /io . 


where I r is the intensity of light reflected from all three boundaries, and the system is 
arranged as in Figure 4. 


light as 


Approximating as before, we may write the reflected intensity of normally incident 
I r = IoB| + I«(l - K,) ! E 2 e' 2 “ ld > 

4 I„< 1 - R,) 2 ( 1 - R 2 ) 2 R,e- 2 “^- 2 “‘ d * , (8) 


and by neglecting terms that contain RjRj we have 
R „ R, 4 R^- 20 ' 111 4 


Now, if the system is arranged as in Figure 5, which will be called the primed 
system, a similar expression may be derived for R’, in which 


R' = y/io ; 


that is, 
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Figure 5. Film-substrate system in primed mode. 

Noting. that for normal incidence the Fresnel formulas for reflection at a boundary, 
as derived in Appendix B, reduce to 


Ri = (n - l) 2 /( n + l) 2 

R 2 = (n - 2 /(n + 2 , (11) 

R 3 = (n t - l) 2 /(n! + l) 2 
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where 1, n, and nj are the refractive indices of air, film, and substrate, respectively, and 
solving equations (9), (10), and (1 1) simultaneously for n yields 


p 2 (l - z 4 ) + p 2 (R - R 3 z 2 p 2 ) - z 2 ( R* - R 3 ) 
n2 ‘ 2n p 2 (l - z 4 ) - p 2 (R - Rjz 2 ? 2 ) + z 2 (R* - R 3 ) + 1 = ° ’ ° 2) 


where 


_0! l d l a “M2 

z = e 1 and p = e 


For a thick substrate, R 3 terms may be neglected and equation (12) becomes 


? „ 1 - 
n 2 - 2n 


•f R - R 1 z 2 /p 2 
1 - z 4 - R + R r zVp 2 


+ 1 = 0 


(13) 


Making some further approximations, we note that for small ajdj , equation (9) 
may be written as 


R = R t + R 2 + R 3 - ( aa^Rj + 2Q' 1 d 1 R 3 + 2oi2d2R 3 ) , 


(14) 


which may be further written, for thin (d < 100 nm) and not too strongly absorbing films 
on a thick substrate, as 


R — Rj + R 2 + R 3 


(15) 


Thus, putting (1 - Ri )(1 - R 2 ) = 1 - R x - R 2 and using equation (15), we have 


(1 - R t )(l - R 2 ) * 1 - R + R 3 


(16) 


Substituting equation (16) into equation (7) yields 

7 = -Ml ~ !i( 1 ~ R s) 

I 2 ( 1 - R + R 3 ) * 


(17) 
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Now, we may calculate n and k for the film from transmission and reflectivity 
data for monochromatic light of a known wavelength X 0 if the film thickness di and the 
substrate refractive index n] are known. We proceed as follows: 

a. Measure the transmission of the bare substrate ; this yields I 2 . 

b. Calculate R 3 using ni in equation ( 1 1 ). 

c. Calculate p using equation (6). 

d. Measure the total reflectivity of the film-substrate system in the primed and 

unprimed modes as shown in Figures 4 and 5 ; this yields R and R\ 

e. Measure the transmission of the film-substrate system; this yields L . 

f. Use the now-known values of I 2 , I 2 , R 3 , and R in equation (17) to 
obtain z. 

g. Use z, p, R, and R’ values in equation (12) to obtain a value for n. 

h. Use the known dj and the calculated z values in the definition of z to 

obtain a value for a, . 

i. Use aj =47rnK/\ 0 to obtain k. 

In this manner, the Tekucheva method yields quick, rough values for n and k which may 
be used as first approximations in more sophisticated techniques. 


B. Through Ellipsometry 

The indices of refraction for the cermet films under study were measured on a 
Gaertner Model LI 19 ellipsometer. This instrument and its peripheral equipment include 
the following: an external-filtered mercury light source emitting 546.1 nm monochromatic 
light; the appropriate collimators and shutters; two Nicol prisms, one rotatable about the 
source beam axis and one rotatable about the source beam axis and a vertical axis; a 
quarter-wave plate, rotatable about the source beam axis and situated between the two 
prisms; and a photomultiplier tube with associated power amplifier and digital voltmeter. 
The photomultiplier tube measures the intensity of the light after it leaves the source; 
passes through the collimator, first Nicol prism and the quarter-wave plate; bounces off the 
specimen; and passes through the second Nicol prism. The specimen is rotatable about the 
same vertical axis as is the second prism. All rotational orientations can be read directly to 
within 0.01 deg. With this equipment, changes in the polarization of the reference beam 
upon reflection from the specimen can be measured. 

The classical treatment concerning the change in the state of polarization of light 
upon reflection from a bare surface or a surface with a film on it is the old and very elegant 
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approach of Drude [42] , whose theory and equations are essentially exact. By using 
Fosterling’s theory of reflection and refraction by a plane-parallel homogeneous thin film, a 
derivation of Drude’s equations is given in Appendix B. Drude’s equations cannot be solved 
in closed form, necessitating trial-and-error and iteration methods. Thus, before the 
availability of electronic computers, the routine application of the exact equations was 
practically impossible. However, computer programs to solve the Drude equations now 
exist, and, in this study, the National Bureau of Standards program of McCrackin [36] is 
used. , 

The McCrackin program was rewritten in Fortran 4H as modified for the Xerox 
Data Systems (XDS) Sigma Five computer used in this study, and, after extensive 
debugging efforts, the program worked well in evaluating ellipsometric measurements of 
the cermet films under study. The exact procedure used in this stage of the study is 
outlined in Appendix A. 


IV. RESULTS OF OPTICAL EXAMINATION 


In this section, the results of the examination of cermet films of various thicknesses 
sputtered from targets with different lead contents are presented. 


A. Microprobe Results 

The AMR microprobe attachment gave values of 58 ± 6, 45 ± 5, and 30 ± 3 percent 
lead by volume for films from the 50-, 25-, and 10-percent lead targets, respectively. The 
data for these films are shown in Figures 6 through 9. The peak height of the PbM (1, 2) 
doublet [43] was taken to be directly proportional to the percent lead by volume 
contained in each film. The peak height of the pure lead film in Figure 6 was used as a 
standard. The errors arise mostly from the simple ratio technique used to obtain 
quantitative values, which is not strictly accurate for the complex metal-dielectric matrix at 
hand. The PbLB (1,2) and MB lines are distorted by the proximity of the barium LB1, 
LB2 (with satellites), LB3, LB4, and LB6 lines, and the platinum M3-01 line (arising from 
the sample holder), and may be disregarded. It might be noted here that the dielectric part 
of the films consisted of the oxides of zinc, aluminum, silicon, boron, and traces of barium. 
These compounds are contained in the Coming 7575 and 7572 lead glasses. 3 


B. Electron Microscope Results 

1. Direct Transmission Studies. Examination of the films deposited onto 
Formvar-covered grids revealed that the films from each sputtering target had their own 
characteristic structure, which may be seen in the following sections. 


3. Brent W. Springford, Private Communication, Corning Glass Works, Corning, New York, 
1970. 
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Figure 8. Spectral pattern of film sputtered from 25-percent pure lead target. 
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40 



1 2 0 ANGLE (deg) 

Figure 9. Spectral pattern of film sputtered from 50-percent pure lead target. 

a. Low-Lead Content Films (30 Percent by Volume). Films deposited 
from the 10-percent pure lead target had a structure consisting of small spheres of lead 
embedded in dielectric. The spheres were typically 10 to 30 nm in diameter and had their 
centers separated by about 40 nm. The lead spheres occupied about 15 percent of the Film 
volume. This structure may be seen in Figure 10. The film was stable under the microscope 
electron beam. 
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b. Intermediate-Lead Content Films (45 Percent by Volume). As the 
metallic content increased, the lead balls seemed to grow in size, typically to about 400 to 
800 nm in diameter, although a number of balls less than 100 nm in diameter may be seen 
(Fig. 11). In other photographs, some of the smaller balls are in the process of combining 
with some of the larger balls, and some of the larger balls are intergrowing (Figs. 1 1 and 
12). A higher magnification of one of the larger balls (Fig. 13) shows an outer edge 
consisting of a number of smaller balls the size of those observed in the low-lead content 
films. The larger balls are usually separated by a center-to-center distance on the order of 
1300 nm. 
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Figure 10. A 30-percent lead film on 
Formvar-covered grid (260 000 x). 


Figure 1 1 . A 45-percent lead film on 
Formvar-covered grid (31 500 x). 



Figure 12. A 45-percent lead film on 
Formvar-covered grid (76 000 x). 


Figure 13. A 45-percent lead film on 
Formvar-covered grid (260 000 x). 
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Since the films are only 40 to 50 nm thick and are only 45-percent 
lead, insufficient metal exists in the films to form spheres of the sizes and distributions 
observed. A more likely shape for the metal particles is some type of disk. 

Let us assume that the particles are oblate spheroids and that they 
contain all the lead in the film (45 percent by volume). This shape is not unreasonable and 
does not conflict with the micrographs. 

Furthermore, the high contrast between the spheroids and the 
remaining film in Figure 1 1 leads one to believe that almost all the lead is contained in the 
particles. Let us further assume that the rotated ellipses forming the oblate spheroids have 
the same eccentricity. A rough volume calculation based on the above assumptions implies 
that such spheroids have eccentricities of about 0.96. 

The above model infers that the larger spheroids have a maximum 
thickness of 200 nm. They should therefore be jutting out of the surface as manifestations 
of the underlying metal structure. 

It should be noted here that it is possible that Figures 1 1 through 1 3 
do not represent the totally undisturbed metallic structure. Some lead migration, as a result 
of film heating in the electron beam, was observed; thus, the original film structure may be 
different from that recorded. Such migratory effects have been noted in pure lead films 
[38, 39]. The effect is not surprising since under certain matrix conditions the electron 
beam can heat the sample to temperatures as high as 1000°C [40]. It is difficult to 
estimate the degree of the migration since the effect occurred very rapidly, usually while 
the microscope was being focused. 

It might be further noted that no migration occurred in the high- or 
low-lead content films except immediately before film breakdown. Apparently, in the 
highly metallic films, conduction to prevent extreme heating existed, and in the slightly 
metallic films the relatively small amount of lead was held in place by the heat-resistive 
dielectric. 

An indication of the accuracy of Figures 1 1 through 1 3 is obtained 
in the replica study presented later in this report. 

c. High-Lead Content Films (58 Percent by Volume). In the extreme 
case, the metal particles remain about the same size as in the intermediate case, but more 
numerous and closely packed; they touch or nearly touch and form an almost continuous 
network (Fig. 14). However, there is not nearly as much contrast between the metallic and 
dielectric areas as there was in the intermediate case, which indicates that lead atoms are 
diffused throughout the dielectric. 

In some areas (Fig. 15) the metal particles have combined 
completely, forming a metallic medium containing dielectric particles that are 
approximately 200 nm across. The lack of contrast in these areas indicates a diffused lead 
content. 
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Figure 14. A 58-percent lead film on Figure 15. A 58-percent lead film on 

F ormvar-covered grid (76 000 x). Formvar-covered grid (260 000 x). 


2. Replica Studies. Replicas were made as follows. Impressions were made 

with replica tape of the various films. The impressions were coated with a sputtered-carbon 
film. This carbon film was shadowed with chromium, and the replicating tape was dissolved 
away in acetone. This is a standard technique [41 ] . 

a. Low-Lead Content Films (30 Percent by Volume). Films that were 

20, 40, 60, 80, 100, and 140 nm thick were examined by the replicating technique. No 
surface trend with the thickness was seen. All films seemed to be continuous and exhibited 
pits and mounds. It is possible that some of the features of the thinner films are caused by 
substrate effects. 


Since the limit of replica resolution is about 5 nm, one might expect 
difficulty in spotting surface features caused by the underlying 20-nm-diameter lead 
spheres. However, some surface details that could be so caused can be seen in Figure 16. 



Figure 16. Replica of 20-nm, 30-percent 
lead film (79 300 x). 


b . I ntermediate-Lead 

Content Films (45 Percent by Volume). For 
this case, films with thicknesses of 20, 40, 
60, 80, and 100 nm were examined by the 
replicating technique. Again no trend with 
thickness was established; all films seemed 
to be continuous. However, several features 
can be seen that could be surface manifesta- 
tions of the lead spheroids observed in 
Figures 11 through 13. The surface features 
in Figures 17 and 18 correspond especially 
well in size and distribution to the sphe- 
roids. This is evidenced because the sphe- 
roids are not totally a product of heating in 
the electron beam. The shape of the fea- 
tures implies that the oblate spheroid model 
is correct in some instances, at least. In 
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Figure 17. Replica of 60-nm, 45-percent Figure 18. Replica of 60-nm, 45-percent 
lead film (31 000 x). lead film (63 000 x). 

these instances, most of the lead in the film is contained in the spheroids. It is possible that 
in areas which do not exhibit such sharp surface features the lead is simply more dispersed 
through the dielectric. 

c. High-Lead Content Films (58 Percent by Volume). Films that were 
20, 40, 60, and 80 nm thick were examined by the replicating technique. A trend with 
thickness appeared and can be seen in Figures 19 through 22. By assuming that most of the 
surface features are a manifestation of the underlying metal network, the plates imply that 
the network becomes more defined and connected as film thickness increases. Metal 
particles seem to form and join together as the film thickens. Definite pits and protrusions 
were seen (Fig. 23) which matched well with the sizes and distributions of the metal 
particles observed in Figure 1 4. 



Figure 19. Replica of 20-nm, 58-percent Figure 20. Replica of 40-nm, 58-percent 
lead film (3 1 300 x). lead film (33 600 x). 
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Figure 21 . Replica of 60-nm, 58-percent Figure 22. Replica of 80-nm, 58-percent 
lead film (30 500 x). lead film ( 1 76 000 x). 



Figure 23. Replica of 60-nm, 58-percent 
lead film (100 000 x). 


C. Interpretation of Electron 
Microscopy Observations 

1 . Internal Film Structure. 
Direct transmission electron microscopic 
examination of lead-dielectric cermet films 
indicates that increases in the lead content 
of the films from 30 to 45 to 58 percent 
vary the film structure drastically. In the 
30-percent lead films the structure consists 
of 10- to 30-nm-diameter lead spheres 
embedded in a slightly homogeneous die- 
lectric medium. As the lead content is 
increased to 45 percent, the spheres grow 
to well-defined oblate spheroids which are 
400 to 800 nm in diameter. Further addi- 
tion of lead (to 58 percent) results in many 
additional spheroids in the 400- to 
800-nm-size range that join to form a con- 
necting network as the film thickens. 


2. Film Surface Structure. Replicas of film surfaces reveal mounds and pits 
that are roughly the size of the underlying metal particles. These features can be 
interpreted as manifestations of the film’s metal structure, which juts out of the film 
plane in some areas and depresses the film surface in other areas. Thus, there appears to 
be a gratifying correspondence between the replica and direct electron microscopy 
studies. 
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3. Effects of Film Thickness Variations. No outstanding thickness trend 
appeared in the range observed except for the case of the 58-percent lead films, which 
seemed to show an increase in the definition of their metallic structures as film thickness 
increased. All films were continuous at 20 nm. 


V. RESULTS OF MEASUREMENT OF OPTICAL PROPERTIES 

This section presents the results of measurement of optical properties of thin 
cermet films. 


A. Refractive Indices from Reflection and Transmission 
Measurements 

Since the noise level in the Cary Model 14 recording spectrophotometer was higher 
than the differences in reflection and transmission that were measured in the less absorbent 
films, physically meaningful refractive indices were obtained only for the 58-percent lead 
films. These values are presented in Table 1 and in Figures 24 and 25. It is difficult to 
estimate the error in these indices; it is merely known that it is large. It should be noted 
that the usefulness of these rough values was in their specification of the initial iteration 
range to be used in the McCrackin computer program; in this sense, they were extremely 
valuable. It is interesting that the rough values for n indicate the same thickness 
dependence that is obtained from the final ellipsometric values. 

TABLE 1. REFRACTIVE INDICES FROM REFLECTION 
AND TRANSMISSION MEASUREMENTS 


Lead Content of Film 
(Percent by Volume) 

Film 

Thickness (nm) 

n'(Real Part 
of Index) 

Kappa («') 

58 ± 6 

80 ± 4 

3.9 ± ? 

0.60 ± ? 

58 ± 6 

60 ± 3 

4.3 ± ? 

0.30+ ? 

58 ± 6 

40 ± 2 

5.0 ± ? 

0.50 ± ? 

58 ± 6 

20+ 1 

2.5 ± ? 

0.60 ± ? 


B. Refractive Indices from Ellipsometry 

Successive iterations of decreasing range and increment values on the XDS Sigma 
Five and 3 hours of computer time yielded the refractive indices shown in Table 2 and 
Figures 26 through 31. The data used are shown in Table 3. Figures 26 through 28 indicate 
a thickness dependence for the real part of the refractive index. No analysis of this 
dependence will be attempted in this study. 
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Figure 25. Approximate values of k' from reflectivity and transmission 
measurements for 58-percent lead films. 
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TABLE 2. REFRACTIVE INDICES FROM ELLIPSOMETRY 


Lead Content of Film 
(Percent by Volume) 

Film 

Thickness (nm) 

n'(Real Part 
of Index) 

Kappa (k') 

58 ± 6 

80 ± 4 

1.94 ± 0.06 

0.35 

± 0.02 

58 ± 6 

60 ± 3 

2.75 + 0.08 

0.51 

± 0.06 

58 ±6 _ 

40+2 

2.69 ± 0.20 

0.79 

+ 0.08 

58 ± 6 

20+ 1 

1.81 ± 0.08 

0.99 

+ 0.05 

45 ± 5 

80 + 4 ^ 

2.15 + 0.02 

0.10 

+ 0.03 

45 ± 5 

60 ± 3 

2.31 ± 0.01 

0.08 

± 0.02 

45 ± 5 

40 ± 2 

2.30 ± 0.02 

0.06 

± 0.02 

45 + 5 

20+ 1 

2.13 ± 0.03 

0.05 

± 0.02 

30 ± 3 

80 ± 4 

2.05 ± 0.02 

0.17 

± 0.02 

30 ± 3 

60 ± 3 

2.16 ± 0.02 

0.05 

+ 0.02 

30 ± 3 

■ . . 40 ± 2 

2.13 ± 0.02 

0.04 

± 0.01 

30+3 

20 + 1 

1.98 ±0.03 

0.020 

± 0.004 


Note. n complex ~ n ' n K * 
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Figure 27. n ' versus film thickness for 45-percent lead films. 
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Figure 28. n' versus film thickness for 30-percent lead films. 



FILM THICKNESS (nm) 

Figure 29. k ' versus film thickness for 58-percent lead films. 
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TABLE 3. ELLIPSOMETER READINGS 


1. Typical Bare Substrate (Kodak B351 Cover Glass) 


Angle of Incidence: 55 deg 

a. Quarter Wave Plate : +45 deg 

Polarizer (deg) 

Analyzer (deg) 

48.85 

178.49 

139.38 

2.09 

b. Quarter Wave Plate: -45 deg 

Polarizer (deg) 

Analyzer (deg) 

47.72 

178.09 

137.33 

2.18 

2. 80-nm, 58-Percent Lead Film on Kodak B351 Cover Glass 

Angle of Incidence: 55 deg 


a. Quarter Wave Plate : +45 deg 

Polarizer (deg) 

Analyzer (deg) 

155.09 

21.44 

65.02 

161.60 

b. Quarter Wave Plate : -45 deg 

Polarizer (deg) 

Analyzer (deg) 

30.45 

160.28 

120.39 

20.53 
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TABLE 3. (Continued) 


3. 60-nm, 58-Percent Lead Film on Kodak B351 Cover Glass 


Angle of Incidence: 55 deg 


a. Quarter Wave Plate: +45 deg 


Polarizer (deg) 

Analyzer (deg) 

148.44 

32.45 

58.20 

152.67 * 

b. Quarter Wave Plate: -45 deg 


Polarizer (deg) 

Analyzer (deg) 


36.65 151.89 

126.78 31.59 


40-nm, 48-Percent Lead Film on Kodak B351 Cover Glass 


Angle of Incidence: 55 deg 


a. Quarter Wave Plate : +45 deg 


Polarizer (deg) 

Analyzer (deg) 

59.36 

150.40 

149.48 

34.81 

b. Quarter Wave Plate: -45 deg 


Polarizer (deg) 

Analyzer (deg) 


125.72 

35.61 


33.82 

149.46 










TABLE 3. (Continued) 


5. 20-nm, 58-Percent Lead Film on Kodak B351 Cover Glass 


Angle of Incidence : 55 deg 


a. Quarter Wave Plate : +45 deg 


Polarizer (deg) 0 

Analyzer (deg) 

151.70 

27.39 

61.54 

156.93 

b. Quarter Wave Plate : -45 deg 


Polarizer (deg) 

Analyzer (deg) 

33.22 

156.09 

123.20 

26.51 

6. 80-nm, 45-Percent Lead Film on Kodak B35 1 Cover Glass 

Angle of Incidence: 55 deg 


a. Quarter Wave Plate: +45 deg 


Polarizer (deg) 

Analyzer (deg) 

51.37 

155.28 

141.53 

17.10 

b. Quarter Wave Plate: -45 deg 


Polarizer (deg) 

Analyzer (deg) 

123.65 


33.08 









TABLE 3. (Continued) 


7. 60-nm, 45-Percent Lead Film on Kodak B351 Cover Glass 


Angle of Incidence: 55 deg 

a. Quarter Wave Plate: +45 deg 

Polarizer (deg) 

Analyzer (deg) 

139.32 

24.48 

49.20 

159.23 

b. Quarter Wave Plate: -45 deg 

Polarizer (deg) 

Analyzer (deg) 

45.85 

158.48 

135.60 

25.00 

8. 40-nm, 45-Percent Lead Film on Kodak B35 1 Cover Glass 


Angle of Incidence: 55 deg 

a. Quarter Wave Plate : +45 deg 

Polarizer (deg) 

Analyzer (deg) 

130.70 

22.28 

39.88 

160.66 

b. Quarter Wave Plate: -45 deg 

Polarizer (deg) 

Analyzer (deg) 

54.30 

160.77 

144.64 

23.26 
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9. 20-nm, 45-Percent Lead Film on Kodak B351 Cover Glass 

Angle of Incidence: 55 deg 


a. Quarter Wave Plate: +45 deg 


Polarizer (deg) 

Analyzer (deg) 

31.83 

. 167.09 

120.20 

13.22 

b. Quarter Wave Plate: -45 deg 


Polarizer (deg) 

Analyzer (deg) 

152.47 


64.31 


10. 80-nm, 30-Percent Lead Film on Kodak B35 1 Cover Glass 


Angle of Incidence: 55 deg 

a. Quarter Wave Plate: +45 deg 

Polarizer (deg) Analyzer (deg) 


57.98 163.21 

148.14 20.26 

b. Quarter Wave Plate: -45 deg 

Polarizer (deg) Analyzer (deg) 


126.86 

37.93 


19.73 

162.76 









TABLE 3. (Concluded) 


11. 60-nm, 30-Percent Lead Film on Kodak B351 Cover Glass 

Angle of Incidence: 55 deg • 

Quarter Wave Plate: -45 deg 


Polarizer (deg) 

Analyzer (deg) 

46.85 

161.05 

136.64 

22.14 

1 2. 40-nm, 30-Percent Lead Film on Kodak B35 1 Cover Glass 

Angle of Incidence: 55 deg 


Quarter Wave Plate: -45 deg 


Polarizer (deg) 

Analyzer (deg) 

144.04 

20.20 

54.55 

163.46 

13. 20-nm, 30-Percent Lead Film on Kodak B351 Cover Glass 

Angle of Incidence : 55 deg 


Quarter Wave Plate : 45 deg 


Polarizer (deg) 

Analyzer (deg) 

65.57 

169.45 

154.94 

11.79 














Figure 31 . k ' versus film thickness for 30-percent lead films. 
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VI. REFRACTIVE INDICES FROM THE MAXWELL GARNETT THEORY 


The theory of Maxwell Garnett [ 27 ] describes the optical properties of a medium 
containing minute metal particles. Maxwell Garnett assumed that the metallic particles 
were spherical and small with respect to the wavelength of light. His results agree with the 
more general theory of Mie [29] for light absorption and scattering by a sphere of arbitrary 
size in the limit of small spheres far apart. 

The Maxwell Garnett theory has received considerable elaboration in the hands of 
Schopper [30] , who assumed nonspherical particles, but it was later shown by Doremus 
[32, 33] that it was not ‘necessary to assume that the particles were nonspherical or to 
make other modifications of the theory. Doremus [32, 33] obtained good agreement with 
experimental data for thin island gold films by use of the Maxwell Garnett theory. 

The theory of Maxwell Garnett will be used in this study to derive theoretical 
values which will then be compared to the experimental values already obtained. The 
method of Male [44] will be used to obtain optical constants from the Maxwell Garnett 
theory. 


In the Maxwell Garnett theory [27] , the effective complex refractive index for the 
film is related to the complex refractive index of the metallic particles contained in the film 
by: 


n ' 2 


u 

+ 2v t 


n 


u 


rr + 2v 


7 


(18) 


where 


n' = n'(l - iO = effective complex refractive 
index for this film 

M = volume factor = fraction of volume taken up by 
the metal particles 

n = ft(l - k) = complex refractive index for the 
metal particles 

v = refractive index of the dielectric medium. 

Equation (18) is the basic equation of the Maxwell Garnett theory and is derived in the 
following subsection. 


A. The Maxwell Garnett Theory 

Maxwell Garnett [27] dealt with the problem of light passing through a dielectric 
medium that contained many small metallic spheres in a volume with linear dimensions on 
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the order of wavelength. In this classical treatment. Maxwell Garnett showed that such an 
assembly is equivalent to a medium with a certain complex refractive index n' = n'(l - 
k' ), and he obtained expressions for fi' in terms of the optical characteristics of the 
metallic spheres and the dielectric medium. His expressions will be derived in this section. 

We shall first need to derive the Lorentz-Lorenz formula. Its derivation follows the 
treatment of Bom and Wolf [45] . 

Let us consider the electromagnetic field in a nonconducting medium in a region 
where the current and charge densities (J and p, respectively) are zero. Using the 
rationalized mks system of units. Maxwell’s general equations are: 

V x H = J + 9D/9t 

V x E = -9B/9t , 

(19) 

V • D = p , 

V • B = 0 , 

where 

E = electric vector 
B = magnetic induction 
J = current density 
D = electric displacement 
H = magnetic vector. 

The last three vectors are needed to describe the effect of the field on material objects. 

We describe the interaction of field and matter through the following relations: 

j = O'E , 

D = e 0 E + P , (20) 

B = (5 + M )p 0 » 
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where 


a = specific conductivity of medium 
P = electric polarization of medium 
M = magnetic polarization of medium 
e 0 = permittivity of free space 
Mo= permeability of free space. 

The terms P and M vanish in a vacuum. 

For the medium under consideration, we have 


o = 0 

and 


j = p = 0 . 

(21) 

Substituting equations (20) into equations (19) and imposing the conditions of 
equation (21) yield 

— 3E — — 

V x B - Mo e o"g£ = Mo(^ x M + 3P/3t) , 

= JMo 

(22) 

where 


j = 3P/3t + V x M , 

(23) 

V x E + 3B/3t = 0 , 

(24) 

_ - V • P o 

v • e = - - , :: 

e 0 e 0 

... •: (25) 


where 


P = - V * P 

= free charge density, 
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and 


V • B = 0 . (26) 

Since the divergence of the curl of any vector is. zero, equation (26) is satisfied by 

B = curl A , (27) 

where A is an arbitrary vector function of position and time, Substituting equation (27) 
into equation (24) and interchanging the curl and time derivation operators yield 

V X (E + 3A/ at) = 0 ; (28) 

and since v X V <t> = 0, where 0 is an arbitrary scalar function, equation (28) is satisfied if 

E = - 3 A /3t - grad<£ , (29) 

Now A and 0 , which are called the magnetic vector potential and the electric 
scalar potential, respectively, must be defined to satisfy the remaining Maxwell equations 
[equations (22) and (25)] . 

Substituting equations (27) and (29) into equations (22) and (25) and noting that 

V x ( V x z) = V ( V • Z) - V 2 Z and ( V • V )Z = V 2 Z , 
where Z is any vector, we obtain 

V 2 A - p 0 € 0 3 2 A/9t 2 - V(V • A + n 0 e 0 d<p/d t) = - M 0 j 5 

(30) 

V 2 0 - jj, 0 e (jd^/dt 2 + 3/3t(v • A + p o e o 30/3t) = - p/e 0 
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If we now impose the Lorentz condition 


V • A + n 0 6(fi<t>/dt = 0 


( 31 ) 


equations (30) reduce to the inhomogeneous wave equations 
V 2 A - n 0 e 0 d 2 A/dt 2 = - Mol ; 

- M o e o 3 2 0/at 2 = - p/e 0 

We should note that since V • v X Z = 0 , equation (22) requires that 


(32) 


- V • 3E/at = V • 7/e 0 


(33) 


and by using equation (25), equation (33) becomes 


dp/d t + V • 7 = 0 


(34) 


Equation (34), the equation of continuity, is satisfied by p and j and may be seen from 
the definitions of p and J [equations (23) and (25)] by applying the vector 
identity V ■ V X Z = 0 . 

The solution of equations (32) may be expressed in terms of the polarization and 
magnetization in the form 

0(r, t) = - - J ->/(l/R)[div» P]dV' , (35) 

47re 0 J 


and 


A (r, t) = /( 1/K) [ curl' M + 9P /at] dVV , (36) 


where the differential operators div ' and curl' are taken with respect to the coordinates (x', 
y ', z') of the variable point of integration r' at which the volume element dV'is situated. 
The square brackets denote retarded time values; i.e. , t is replaced by tA/Mo^o R inside 
each bracket. Further, 
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R = | r - r* | = {(x 


(37) 


- x ') 2 + (y - y') 2 + ( z 


z ’) 2 } 


7 


2 


is the distance from the point r(x, y, z) to the volume element dV' at r'(x y z'). The 
integration is performed over all space. 

Let us verify that equation (35) is a solution for the scalar, potential in equation 
(32). Imagine that the point r is surrounded by a sphere of radius a, centered on r . Divide 
equation (35) into two parts: 


0-01 +02 


(38) 


where 


<p = contribution to the integral from the interior of the sphere 
02 = contribution from the rest of space. 

Since R 0 for every point r'(x', y', z') , there are no singularities in 

the 02 contribution, and 0 2 may be differentiated under the integral sign. 

It should be noted that equation (35) may be written 

0(r, t) = /p(r, t - 'sT^T 0 R)/e dV* . (39) 

In spherical coordinates when considering solutions representing spherical waves, i.e., 
solutions of the form V = V(r, t), we have 


V 2 V = (l/r 2 )a /9 r(r 2 9 V/3 r) 


(40) 


and we may write 


V 2 (p/r) = ( l/r)9 2 p/9r 2 


(41) 


thus showing that <t> 2 satisfies the homogeneous wave equation 


V 2 0 2 “ ^e$ 2 <p 2 /dt 2 = 0 


(42) 
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since 


V v 2 - *,»«*/*» . i 8 , i& t ^ r) 


Po e o 9 2 p(r, t - Vp^r) A 

“T " 0 


( 43 ) 


Now the case of <t> i will be considered, noting that the integrand has a singularity 
at R = 0. It will be assumed that p' is a continuous function of r and t. 


By making the radius of the sphere sufficiently small, we can assure that for 
all r ' inside the sphere 


p( r* , t - i ^T 0 B) - p( r, t) < 8 


(44) 


for any prescribed 8 , where 8 > 0. Thus, as the radius => 0 1 V 2 <fi i approaches more 
closely the electrostatic potential of a homogeneously charged sphere with charge 
density p ; i.e., 

V 2 0i = - p(r, t)/e 0 (45) 


for sufficiently small a. 

Equation (45) is called Poisson’s equation; it arises directly from the definition 
of 0 and the electrostatic field equations. 

Also, as a -> 0, 0j => 0, since if a is sufficiently small, we may write 

• • •• IT 



i lh t ) 

47T60 


/ ( 1/R) dV» 

R<a 


_£ — r 2 it r 2 f a r gin q (j r 

2ire 0 i i i 



Therefore, it follows that as a -*■ 0, 


0 as a — 0 


(46) 


V 2 0 " Mo e o0 = V 2 (^>i + <pj - Mo e o (0i + <P 2 ) p(r, t)./e 0 . (47) 
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Thus equation (35) is a solution for the scalar potential in the inhomogeneous wave 
equation (32). In an exactly analogous manner, it can be seen that each Cartesian 
coordinate of equation (36) is a solution of the corresponding inhomogeneous wave 
equa tion, noting that in the process we would substitute for p' the appropriate component 
of V Mo To j. Thus, equation (36) satisfies the inhomogeneous wave equation (32) for the 
vector potential. 

Note that since p" and j satisfy the continuity equation (34), equations (35) and 
(36) may be seen to satisfy the Lorentz condition of equation (31). 

One may regard A and <t> in equations (34) and (35) as arising from contributions 
from each volume element in space. Note that s/no^o R is the time required for light to 
travel from r' to r . Equations (34) and (35) are therefore called retarded potentials. 

We now use the following identities: 

div’ [ P ] = [ div’ P ] + r . [ p] 

and (48) 

curl* [ M] = [ curl* M ] + R x [ M] , 

•K 

where 


R = r -F , 

to write equations (35) and (36) as 

* =- ^/{i dlv ’ [5i - • ^i} dv V 

A = ■— fU curl’lM] - R x [M] + [ -p[ dv ’ • 

4w J | R R R t 

Equation (46) will now be rewritten fey modifying the vector identities 

div(P/R) = (1/R) div P + P • grad(l/R) 
and 

curl (M/R) = (1/R) curl M - M X grad(l/R) . 

According to the divergence theorem, 


(49) 


(50) 
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/ 


/ c • n dS = J\v • C dV , (51) 

where C is any vector field C(r) defined within a volume V surrounded by the closed 
surface S . 

Equation (51) may be used to write equations (50) as 

/( 1/R)(n • P) dS = f { ( l/R) divP + P • grad (l/R)} dV 


and (52) 

J ( l/R) (n x M)dS = f { (l/R) curl M - M xgrad(l/R)} dV, 

where we have integrated over an arbitrary finite domain, letting h be the unit outward 
normal to the surface. 

Assume now that the region of space where P and M differ from zero (the 
substance) remains within a finite closed surface. If the integrals in equations (49) are 
extended over the volume inside this surface, the surface integrals in equations (52) vanish, 
and equations (49) can be written as 

<P = /{[P] • grad’ ( l/R) + R • [P] j.dV’ 

01 ' ' (53) 

A = /j [M] x grad’ ( l/R) - R x [M] + [P]/R jdV’ . 

Let us now consider the substance itself. We shall regard matter as being composed 
of interacting particles embedded in the vacuum. As long as the region over which our 
averages are taken is large compared with the particle dimensions, the electromagnetic 
properties of each particle may be simply described by electric and magnetic dipoles. 

We shall now consider a small sphere centered on a particular molecule, with a 
radius that is large compared to the linear dimensions of the molecule. Neglecting 
molecular structure outside the sphere we assume that the polarization P is constant. We 
shall assume that the molecules inside the sphere cause no resulting field at the central 
molecule. We therefore have a molecule situated in a spherical region, inside of which there 
is vacuum and outside of which there is a homogeneously polarized medium. 

Now, 4> must be determined for this case, namely <t> for the free charges on the 
spherical discontinuity surface where P changes from zero (interior) to a constant exterior 
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value. Consider the potential 0 of the complementary configuration, that of a 
homogeneously polarized sphere surrounded by vacuum. The superposition of these two 
configurations is a homogeneously polarized surface with no boundary. Therefore, the 
potential caused by the boundary is zero and 


0 + 0=0 


(54) 


The term 0 is just the scalar potential stated in equation (53), but with P being a 
constant; thus, 


0 = - 0 = - ~~ P * f grad* ( l/R) dV* . (55) 

4tT€q j 

Now since R = {(x - x') 2 + (y - y') 2 + (z - z') 2 } ^ , grad' can be replaced in equation (55) 
by — grad and equation (55) can be written as 


0 = — 0 = 


1 — 


47re n 


P * grad J ( l/R) dV* 


- P * grad0 o 


(56) 


where 


00 * ' 4^ ( (l/E)dV ’ • (57) 

The term 0 O is interpreted as the potential of a uniformly charged sphere of 
charge density p = - 1 . Then by Poisson’s equation. 


V 2 0 o = l/e„ 


(58) 


Now, by using equation (56), 

-00/ax = ^-(P 30 o /9x + P 90 o /9y + P 90 o /9z) 
ctA x y z 


= P 9 2 0 o /9x 2 + P 9 2 0 o /9x9y + P 9 2 0 o /9x9z ■ ' (59) 

x y z 
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with similar expressions for 30/3y and -30/d z . Now, by symmetry it is known that at 
the center of the sphere 


3 2 0 o /9x9y = 9 2 0 o /3x3z = d^Q/dydz = 0 


'(60) 


and 


a V 9 * 2 = aV3y 2 = d 2 h/^ 2 = » (61) 

je 0 

where equation (58) has been used in equation (61 ). 

By using equation (59) and similar expressions and noting equations (60) and (61), 

-vV> = f p a 2 0 o /ax 2 + f p a 2 0 o /ay 2 + £p a 2 0 o / 9z2 = • (62) 

x y z 3e 0 

The total field inside the sphere, E ', which is the effective field acting on the 
central molecule, is obtained by adding the mean field E to equation (62), which yields 

E» = E + . (63) 

J€ o 

It will be assumed that for each molecule the electric dipole moment p established under 
the influence of the field is proportional to the effective field E' ; i.e., 


p = «E’ 


(64) 


where a represents the mean polarizability. Thus, with N molecules per unit volume the 
total electric moment P per unit volume is 


P = Np = NaE’ . (65) 


Assuming that P is proportional to E, that is, . 


P = tjE 
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( 66 ) 



where t? = dielectric susceptibility, we'have, by eliminating E' between equations (63) and 
(65), that f . . 


77E 

Na 


E + 


jE 

3e 0 


or 

Na 

V ■' 

" 3e 0 

Now using equation (20) 

D = e 0 E + P = eE = e 0 E + 7]E , 

where an isotropic medium has been assumed with 


D = eE 


* 


where e == permittivity. Thus, from equation (69) 
e = e 0 + 77 

Substituting equation (68) into equation (71) yields 
2Nc* 

£ ° + 3 

C Ka_ . 

3e 0 

and solving equation (72) for a yields 


3 £ o / 


w°| 

CO 1 
11 

,n l " : ) 


N V 

» e + 2e n / 

N ' 

v n 2 + 2/ 

> 


(67) 


( 68 ) 


(69) 


(70) 


(71) 


(72) 


(73) 
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where the Maxwell relation e/e 0 = n 2 for nonmagnetic substances has been used. 
Equation (73) is called the Lorentz-Lorenz formula. 

Equation (73) will be used to obtain the Maxwell Garnett formulation. Note first 
that by using equations (63), (65), and (73), 

2 

E* = E + ~ i E’ (74) 

+ 2 

Let us now consider light of a wavelength X falling on a metal ball with 
radius a and optical constants n , the real coefficient of refraction, and x , the 
absorption index. The total field inside the sphere is then given by equation (74), subject to 

n = ft(.l - be) • (75) 


If we now consider light falling on a medium consisting of many small metal spheres 
distributed in a vacuum, many to a wavelength, and integrate over a unit volume, this has 
the effect of summing the field contributions of the spheres, and equation (74) becomes 


E f = E + (dirN’aVs) 


- 1 

n 2 + 2 


(76) 


where 


N' = number of metal spheres per unit volume of medium 

4TrNa 3 /3 = percent metal per unit volume of medium. 

Now E' represents the effective field in the new medium. It has been assumed that the 
medium under consideration extends throughout a space of dimensions which in no 
direction are smaller than the order of a wavelength. 

Thus we have a medium with a revised polarization P' , which by definition means 

that 


_ _ p» 

E’ = E + 


3e 


o 


and by comparing equations (76) and (77) we have 
P' =4,N-a%(^i)E- . 


(77) 


(78) 
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By eliminating E' from equations (77) and (78) and assuming that 


P' = V E , 


where t? ' is the dielectric susceptibility for new medium, we have 


rj* E /n 2 + 2 \ — n' E 

— ,(7rr) - E + 


47rN f a 3 e 


3e n 


or 


,» - 


(4*N-a’€,)->(^f) * ~k, 


(79) 


The dielectric constant e ' for the medium is defined as 


e’ = e 0 + V 


(80) 


Substitution of equation (79) into equation (80) yields 


4,rN ' a 3 €„ (frrr) 


47tN’ a 3 

/n 2 - IN 

• 

3 1 

W + 2 / 



Now, if q=4n-N'a 3 /3, 




e o '+ 


3^e 0 


1 - M 


n 2 - 1 
n 2 + 2 
n 2 - 1 
n 2 + 2 


(81) 


(82) 


If the metal spheres are situated in glass of refractive index instead of a vacuum, eo -*• 
e" in equation (82), where e" is the dielectric constant of the glass, yielding 



(83) 
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where 


n 2 = e/e 0 -*■ e/e" . 


Equation (83) may be written as 


e' e" 


e o 



(84) 


Now by noting that 


e"/e 0 = v 2 , 


and defining 


(n 1 ) 2 = [n'(l -k')] 2 =e'/e 0 , 


we have 


= 


(n f ) 


u 2 + 


2 2 

O 2 n - ” 

3/iU 2 o 2 
+ 2ir 

n 2 - u 2 

1 " M n 2 + 2u 2 


(85) 


where n ' is the effective refractive index of the new medium. Note that equation (85) may 
be rewritten as 


n’ 2 - u 2 
n T 2 + 2u 2 


= 


n 2 - u 2 


n * + 2v 


T 


( 86 ) 


\ 
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since from equation (86) 


5 ’ 2 


9 „ 9 n 2 - u 2 

u 

n 2 - v 2 
1 — M 9 2 

n 2 + 2u 2 


= ir + 


Q 2 n 2 - u 2 
2 3MU n 2 + 2u 2 


2 T 
m - v 


1 ' M n 2 + 2u 2 


Equation (86) is the Maxwell Garnett relation which describes the effective 
refractive index, n' , of a dielectric embedded with many metallic spheres per wavelength 
of light in terms of the optical constants of the dielectric and the individual spheres. It 
should be noted that n' depends on p , the relative volume of metal, rather than on the 
individual sphere radii. Thus, the spheres may be of quite various radii, provided that they 
are present in sufficient numbers. 


B. Results Using the Maxwell Garnett Theory 

Now, following the method of Mal£ [ 44 ] , we proceed as follows. Rearranging 
equation (86) yields 


n’ 


' (n z - u ; ) 

1 " M (n 2 + 2u 2 ) 


2 „ 9 (n 2 - u 2 ) 

- u = 2f,u G7T 2 J) 


( 87 ) 


or 


n’ 2 ( 1 - pr) - u 2 = 2u 2 pr , 


( 88 ) 


where, in polar form 


r = 


n 2 - u 2 
n 2 + 2u 2 


= Ae 


10 ! 


( 89 ) 
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Solving for n' 2 in equation (87), we have 



where the definition-of n ' has been noted and rewritten in polar form. 

Now let us define the following points in the complex plane: 

N = -l/2r D = 1/r Q = p . (92) 

The graphical solution for h' = uBe 1 ^ can now be written since 

B = (2QN/QD) 1/z and 0 - — . (93) 

Further, since h' = n' - in'« ' , 


n’ = -uB cos/3 and n f /c T = uB sin 0 


(94) 


Let us now evaluate r for our case, letting v = 1.5, ft = 2.42, and nic = 4.45 [46] , where 
we have used the definitions in equation (75). From equation (89), it is found that 


r ~ + 2 ^ 


Therefore, by definition of n we may write 



( 95 ) 
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Now substituting numerically into equation (95), 

= -16.2 - 21.6 i _ (-16.2 - 21. 6i) (-9.45 + 21. 6i) 
r ~ -9.45 - 21.6 i (-9.45 - 21. 6i) (-9.45 + 21. 6i) 

= 1. 12 -0. 262 i = Ae® 

1 / 

A = {( 1. 12) 2 + ( 0. 262) 2 } /2 and a = tan -1 (-0.262/1.15) 
= 1.15 


Now, by using equation (92), 

N = -l/2r = -0.434 e' ia , 

D = 1/r = 0.87 e' ia , 

Q = m — 0.15, 0.25, and 0.45 for our three cases, as 
estimated from the electron micrographs. 

By using the plot in Figure 32 and proceeding in accordance with the 
Male", the following cases are obtained: 

Case I. 30 Percent Lead Films 

Q = g = 0.15 . 

Now, from Figure 32, it is found that 

QN = 5.78 QD = 7.26 DQhf = 174.2° . 

Therefore, from equation (93), 

B = ( 2QN/QD) ^ =1.26 0 = + - = 177° 6' 

2 

Now by using these results in equation (94), 

n T = -uB cos/3 = 1.89 n'K 1 = uBsin/3 = 0.096 


( 96 ) 


(97) 


(98) 


method of 


51 



IMG. 

AXIS 


0.6 


0.4 -f 



Figure 32. Male’s method for solution of the Maxwell Garnett equations. 
Similar procedure for the other cases yields the following: 

Case II. 45 Percent Lead Films 
Q = M = 0.25 , 


n' = 2.18 , 


and 


n'x ' = 0.19 . 

Case III. 58 Percent Lead Films 
Q = H - 0.45 , 
n' = 2.91 , 

and 

n'x' = 0.49 . 

The above results are listed in Table 4, along with the extreme values obtained 
experimentally. 
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TABLE 4. THEORETICAL AND EXTREME EXPERIMENTAL 
REFRACTIVE INDICES 


Lead Content of Film 
(Percent by Volume) 

n'(Real Part 
of Index) 

n'/c '(Imaginary 
Part of Index) 


Theoretical 

Experimental 

Theoretical 

Experimental 

30 ± 3 

1.89 

1.95 to 2.18 

0.10 

0.03 to 0.41 

45 + 5 

2.18 

2.10 to 2.32 

0.19 

0.06 to 0.30 

58 ± 6 

2.91 

1.73 to 2.83 

0.49 

0.57 to 2.94 


Note. n com pi ex n "in* 


VII. CONCLUSIONS 


As a result of the physical and optical examinations of the thin lead dielectric 
cermet films studied in this report, one may conclude the following. 


A. Internal Film Structure 

Direct transmission electron miscroscopic examination of lead dielectric cermet 
films indicates that increases in the lead content of the films from 30 to 45 to 58 percent 
vary the film structure drastically. In the 30 percent lead films the structure consists of 10- 
to 30-nm-diameter lead spheres embedded in a fairly homogeneous dielectric medium. As 
the lead content is increased to 45 percent, the spheres grow to well-defined oblate 
spheroids which are 400 to 800 nm in diameter. Further addition of lead (to 58 percent) 
results in many additional spheroids in the 400- to 800-nm-size range which join together 
to form a connecting network as the film thickens. 


B. Film Surface Structure 

Replicas of film surfaces reveal mounds and pits that are roughly the size of the 
underlying metal particles. These features can be interpreted as manifestations of the film’s 
metal structure, which juts out of the film plane in some areas and depresses the film 
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surface in other areas. Thus, there appears to be a gratifying correspondence between the 
replica and direct electron microscopy studies. 


C. Film Refractive Indices 

The refractive indices of lead-dielectric cermet films as measured by ellipsometry 
agree acceptably with the theoretical values predicted by the Maxwell Garnett theory for 
those films whose structure fits that of the Maxwell Garnett model; i.e., for those low- and 
intermediate-lead content films whose structure consists of tiny metal balls embedded in a 
dielectric. As expected, the experimental values for the highly metallic network structure 
films diverge from the Maxwell Garnett values. This is especially noticeable for the 
imaginary part of the refractive index. Thus the Maxwell Garnett theory may be applied to 
thin cermet films of suitable structure. 


D. Effects of Film Thickness Variations 

Electron micrographs reveal no outstanding thickness trend in the range observed 
except for the case of the 58-percent lead films, which seem to show an increase in the 
definition of their metallic structures as film thickness increased. All films examined were 
continuous at 20 nm. 

The ellipsometer measurements reveal a definite thickness trend for the real part of 
the refractive index. This may be seen clearly in Figures 26 through 28. No existing 
theories have predicted such a trend. 
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APPENDIX A. USE OF THE ELLIPSOMETER 


1. Theory of the Ellipsometer 

Ellipsometry is the measurement of the change in the state of polarization of light 
on reflection from the surface of a dielectric or metal covered with a film. The ellipsometer 
in its modern form was developed by Rothen [47] in 1944! Motivated by the need of 
measuring thin film thicknesses quickly and accurately, he designed an apparatus, which he 
called an ellipsometer, based upon the fundamental principles established by Drude [37] 
concerning the change in the state of polarization of light upon reflection from a bare 
surface or a surface with a film on it. Rothen’s apparatus followed the general form shown 
in Figure A-l. With the advent of modern electronic computers the application of the 
somewhat cumbersome exact Drude equations became routine, and ellipsometry became a 
powerful technique for measuring the thicknesses and refractive indices of films. An 
excellent review of the historical development of ellipsometry has been written by Rothen 
[35], 


In understanding the use of the ellipsometer in analyzing polarized light, it is 
convenient to utilize the Poincar^ sphere representation of polarized ‘light. This 
representation will now be developed by use of Stokes parameters, and will, in general, 
follow the presentation of McCrackin, et al. [48] in relating the Poincare sphere 
representation of ellipsometry. 



'-’ire A-l. Schematic representation of the halfshade ellipsometer [46] . 
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Consider a plane monochromatic wave traveling along the axis of a coordinate 
system. For elliptically polarized light, the electric vector of the wave is 

E x = a i cos (t + 6*) , 

E y = a 2 cos ( T + 5 2 ) 


and 


(A-l) 


E = 0 
z 


where 


T = <o(t-z/v) 
co = angular frequency 
v = linear velocity of the light 
aj = component amplitude 

The polarization can be described by the phase difference, 5 = 6 2 - 5, , and the 
amplitudes, aj and a 2 , of the components. In a stationary plane whose normal is parallel 
to the z axis, the locus of the electric vector’s end is the ellipse shown in Figure A-2. 

The ellipse may also be described in relation to coordinates x' and y' along the 
axis of the ellipse. Thus, the semiaxes, a and b , of the ellipse and the 
inclination 0 of these coordinates also describe the polarization of the light. 

Now, the auxiliary angles a and x are introduced and defined by 


tan a = a 2 /aj 


(A-2) 


and 


tanX = ±b/a 


(A-3) 


The angles 0 and x are called the azimuth and ellipticity of the light, respectively; and 
representation may be made either in terms of 0, x and a! + a 2 ; aj , a 2 , and 6 ;or a, 6, 
and a [ + a 2 . 
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Figure A-2. Locus of the electric vector for elliptically polarized light 
in a plane normal to the direction of propagation. 

Still another representation of the state of polarization, and one which leads quite 
naturally to the Poincare sphere, is by parameters which all have the same physical 
dimensions. These are called Stokes parameters and can be defined as 

So = a i 2 + a 2 2 , 

s i = a i 2 - a 2 2 , 

5 2 = 2aja 2 cos 6 , 

and (A-4) 

5 3 = 2a 1 a 2 sin6 

Note that S 0 2 = Si 2 + S2 2 + S3. 2 , so tha't the Stokes parameters are not all independent. 
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The Stokes parameters are also given by 
= S 0 cos 2X cos 20 , 

5 2 = S 0 cos 2X sin 20 , 

and 

5 3 = S 0 sin 2X 

The proof of equations (A-5) follows [45] . 

Returning to equations (A-l), they are rewritten as follows: 

E /a t = cost cos Si - sinr sin 
x 

and 

E^/a 2 = cos r cos 6 2 - sin r sin 6 2 

Hence, 

(E^/aj) sin 6 2 - (E^/a 2 ) sin Si = cos r sin (6 2 - 6 t ) 
and 

(E /aj) cos 6 2 - (E /a 2 ) cos 6i = sinr sin (6 2 - 5 X ) 
Squaring and adding yields 

E E 

(E /a^ 2 + (E /a 2 ) 2 - 2 X — cos 5 = sin 2 5 , 

' x 1 y a i a 2 


(A-5) 


(A-6) 


(A-7) 


(A-8) 
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where 


5 - 6 2 - 


Equation (A-8) is of the form 

Ax 2 + Bxy + Cy 2 + Dx + Ey +, F = 0 (A-9) 

where A, B, C, D, E, and F are constants, which is the general equation of second degree. 
If B 2 - 4AC < 0 , then equation (A-9) represents an ellipse [49] . 

Comparison of equation (A-9) with equation (A-8) yields 

B 2 - 4 AC = (4 cos 2 6)/( ai 2 a 2 2 ) - 4/( ai 2 a 2 2 ) = 4(cos 2 6 - l)/( ai 2 a 2 2 ) 

= -4 sin 2 6/( ai 2 a 2 2 )<0 . (A-10) 


Therefore, equation (A-8) represents an ellipse, and equations (A-l) represent elliptically 
polarized light. The ellipse is inscribed into a rectangle whose sides are parallel to the 
coordinate axes and whose lengths are 2a! and 2a 2 . The ellipse touches the sides at (taj , 
±a 2 cos 5) and (±a! cos 5, ±a 2 ) and has its axis along the coordinates x' and y' (Fig. 5). 
Then, using the standard rotation matrix we have 

E = E cos 0 + E sin 0 
x f x r y r 


and (A-l 1) 

. . ' E = -E sin 0 +. E cos 0 
y* x r y 

Note that 2a and 2b (a < b, Fig. A-2) are the lengths, of the axes of the ellipse, therefore 
E f = a cos ( t + 5 0 ) 

and (A- 12) 

E = ±b sin (r + 6 0 ) > , s - ■ 

where the two signs describe the two possible senses in which the end point of the electrical 
vector may describe the ellipse. 
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Comparison of equations (A-l 1 ) and (A-1 2) and use of equations (A-6) gives 
a (cos t cos 6 0 - sin r sin S 0 ) = a,(cosr cos - sinT sin 6j) cos 0 

+ a 2 (cos r cos 6 2 - sin t sin 6 2 ) sin 0 


and 


±b(sinTCos6 0 + cosrsin6 0 ) = -a t ( cos r cos 6j - sin t sin 6^ sin 0 

+ a 2 (cosTCos6 2 - sinrsinS 2 ) cos <p . 

(A-l 3) 

Now, equating the coefficients of cos r and sin r yields 

a cos 6 0 = a. i cos cos 0 + a 2 cos 6 2 sin 0 , (A-14) 

a sin 5 0 = a A sin 6 X cos <p + a 2 sin 6 2 sin 0 , (A-l 5) 

±b cos 6 0 = a t sin sin <j> - a 2 sin 6 2 cos 0 , (A-l 6) 

and 


±b sin 6 0 = -a 1 cos 6 * sin 0 + a 2 cos 6 2 cos 0 . (A-l 7) 

Squaring and adding equations (A-14) and (A-l 5) and noting that 5 = 5 2 - 5 1 yields 

a 2 = aj 2 cos 2 0 + a 2 2 sin 2 0 + 2a t a 2 cos 0 sin 0 cos 6 . (A-l 8) 


Similarly, from equations (A-16) and (A-17), 

b 2 = a t 2 sin 2 0 + a 2 2 cos 2 0 - 2a 1 a 2 cos 0 sin 0 cos 5 

Hence, 

a 2 + b 2 = aj 2 + a 2 2 


(A-l 9) 


(A-20) 
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Then multiplying equation (A- 14) by equation (A- 16) and equation (A-15) by equation 
(A-l 7) and adding, one obtains 


±ab = a 1 a 2 (cos sin 8 2 - cos S 2 sin 5 X ) = a^ sin 6 . (A-21) 


Further, dividing equation (A-l 6) by equation (A-l 4) and equation (A-l 7) by equation 
(A-l 5) yields 


±b/a = 


a! sin sin 0 - a 2 sin 6 2 cos 0 
a x cos cos 0 + a 2 cos 6 2 sin 0 


-a. i cos 6 1 sin 0 + a 2 cos 5 2 cos 0 
a t sin cos 0 + a 2 sin 6 2 sin 0 


(A-22) 


Rewriting equation (A-22) yields 

cos 0 sin <p (a 2 - a 2 2 ) = a 1 a 2 (cos 6* cos 8 2 + sin 5* sin 6 2 ) 


x (cos 2 0 - sin 2 <p) 
or 


(a 2 - a 2 ) sin 20 = 2aja 2 cos 6 cos 20 


(A-23) 


Now, using a as defined in equation (A-2), equation (A-23) can be written as 
2a 1 a 9 cos 5 2 tan a 

tan 20 = ~~Y n r ~ = , - cos 6 


a i - H 


1 - tan a 


or 

tan (20) = (tan 2a) cos 6 

Now, from equations (A-20) and (A-21), we have 

2ab 2a^ 9 sin 5 _ 2 tan a . _ 

a 2 + b 2 a 2 + a 2 2 1 + tan 2 a SU1 

= 2 sin a cos a sin 5 = (sin 2 a) sin 5 


(A-24) 


/ 


(A-2 5) 



and using x as defined in equation (A-3), equation (A-25) becomes 


sin 2X = (sin 2 a) sin 5 


(A-26) 


Now, from equations (A-4), (A-25), and (A-26), we have 


S 3 = 2a 1 a 2 sin 5 = (a^ + a 2 2 ) sin 2X = S 0 sin 2X 


(A-27) 


Further, it is noted from equations (A-4) and (A-24) that 


S 2 = 2a t a 2 cos 6 = (a t 2 - a 2 2 ) tan 2 0 = S 1 tan 2 0 


(A-28) 


Now, noting that S 0 2 = Si 2 + S 2 2 + S 3 2 and using equations (A-27) and (A-28), we have 
S 0 2 = S 2 + Sj 2 tan 2 20 + S 0 2 sin 2 2X 
which means that 


, _ Sf(l - Sin^ZX) . s 2 2 2 

1 1 + tan 2 20 0 Y 


or 


S} = S 0 cos 2X cos 20 


(A-29) 


Substituting equation (A-29) into equation (A-28) yields 


S 2 = S 0 cos 2X sin 20 


(A-30) 


However, equations (A-27), (A-29), and (A-30) are just equations (A-5), which we desired 
to derive. We can now proceed with the Poincare^ sphere representation of polarized light. 

It should be noted that equations (A-5) are the transformation equations required 
to go from linear to spherical coordinates. Thus, the state of polarization may be described 
by a point on a sphere of radius S 0 by the spherical coordinates S 0 , 20, and 2x . 
Alternately, the Stokes parameters may be used as Cartesian coordinates to represent the 
polarization. This is shown in Figure A-3. 
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With the assumption that the law of conservation of energy holds and that the light 
intensity is unity, the point can be projected onto a sphere of unit radius. The coordinates 
are then given by 


and 


50 - 

_ 1 - tan 2 a. 

51 = : — 5 — = cos 2X cos 20 

1 1 + tan- 2 a 


2 tan o? cos 5 
1 + tan 2 a 


cos 2X sin 20 


( A-3 1 ) 


where S 0 = a! 2 + a 2 2 = 1 and equations (A-2), (A^4), and (A-5) have been used. 

Thus, light of a given polarization is 
represented by a point on the sphere with 
polar angles 20 and 2x , which represent, 
respectively, twice the azimuth and twice 
the ellipticity of the light. Consider an 
example. For plane polarized light, the 
ellipticity x is zero and S 3 = 0 . Thus, 
plane polarized light is represented by 
points on the equator of the Poincare 
sphere. Again, consider circularly polarized 
light, with an ellipticity of 45 deg. Here, Si 
= S 2 = 0 and 2x = ±ir/2,and circularly 
polarized light is represented by the poles 
of the sphere. 

This representation may be used to 
define elliptically polarized light with 
Figure A-3. Representation respect to any physical axis, say x" and 

of polarized light on the Poincar^ sphere, y" at an angle 0" with the xy axes. The 

ellipticity x is independent of the choice 
of axes, but the azimuth 0 is changed by <p" in going from the unprimed to the double 
prime system. Thus the Si and S 2 axes must be rotated about the S 3 axis by 20”. 

The Poincare sphere is convenient for representing the effects of double refracting 
plates and reflection on the state of polarization of a light beam. For example, let the 
polarized light in Figure A-2 pass through a double refracting plate of relative phase 
retardation 6 . Choose the x and y axes in Figure A-2 along the fast and slow axes of the 
plate. Since the azimuth of the fast axis is zero, consideration of equations (A-3 1 ) shows it 
to be represented by the positive Si axis. The only effect of the plate is to change the 
relative phase of the components by 6 ; thus, ai and a 2 remain constant, and by 
equations (A-2) and (A-31), a and S ( also remain constant. Thus, upon passing through 
the plate, the point representing the polarization will remain on the intersection of the 
sphere and the plane Si = constant . Consideration of equations (A-31) shows that this 
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intersection, as illustrated in Figure A-4, is a circle of radius (2 tan a)/(l + tan 2 a) . The 
incident light is represented by the point P , and the transmitted light is represented by the 
point L : Therefore, the effect of the doubly refracting plate is to turn, points on the 
sphere about the Si axis by an angle 6 . If the plate had been oriented with its fast axis at 
a <p' azimuth with respect to the xy axes, the rotation by 9 would have been about the 
line between the center of the sphere and the point 2<p' on the equator. 


s 3 



Figure A-4. Intersection of the Poincare sphere and a plane normal to the S! axis. 

Now, consider the Poincar^ sphere representation of reflection at a metal surface. 
For this case, the incident light ray is resolved into components in the plane of incidence 
and normal to the plane of incidence (in the surface plane). The amplitude components 
Aq p , Aq s of the incident waves and the corresponding components Rq p , Rq s are 

related by the Fresnel relations derived in Appendix B and may be stated as 
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and 


= tan (0p - 9,) 

0,p tan ( 9 0 + 0 t ) 0,p 


(A-32) 


R = sin ( 9 0 - 9 X ) A 
0,s sin (0 O + dj) 0,s 

The quantities involved are further governed by Snell’s law. 


sinP^ (rifl sin 0 0 ) /n t , (A-33) 

noting that since this case ni is complex ;-thus, d x is complex also and no longer has the 
simple significance of an angle of refraction. Further, since 0 x is complex, so are the 
ratios R() p/Ag s and Rq s^O s » he., characteristic phase changes occur on reflection. 

Thus, linearly polarized light will, in general, become elliptically polarized on reflection at 
the metal surface. Then, if 5 p and 5 S are the phase changes and p p and p s , which are 

the absolute values of the reflection amplitude coefficients, we have 


i6 


i6 


R /A n 
0,P 0,p 


= p e 
P 


and 


R /A n 
0,s 0,s 


= p e 
s 


(A-34) 


Further, if the incident light is linearly polarized in the azimuth <p 0 , we have 


tan <t> Q 


= A /A 
0,s 0,p 


(A-35) 


Thus 


i(6 - 6 ) 

p = R /R = (p /p ) ta.n(j)e P S 
0,p 0, s v ^p r s' r 0 


(A-36) 


or, more generally, 


iA 

p = tan ip e 


(A-37) 
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where \p and A are functions of the optical constants of the surface, the angle of 
incidence, and the thickness and refractive index of any surface films present, and A = 5p 

- $ s , representing the phase difference between the s and p components. Tan \ p is simply 
a measure of the relative absorption of the two components. 

For consideration of the reflection process on the Poincare" sphere, represent the 
light with respect to Cartesian coordinates with the z axis in the direction of propagation 
and the x and y axes in and normal to the plane of incidence. Then tan \p = a s /a p . The 

Stokes parameters before reflection are given by equations (A-3 1 ) and after reflection by 


,, _ 1 - tan 2 o? cot 2 tp 

1 1 + tan 2 a cot 2 ip * 

,, _ 2 tan o? cot ip cos (6 + A) 

2 1 + tan 2 a cot 2 ip ’ 

and (A-38) 

_ 2 tan a cot ip sin (6 + A) 

3 1 .+ tan 2 a cot 2 ip * 


where it has been noted that the process of reflection introduces a phase 
difference A between the s and p components and changes the ratio of their amplitudes 


by tan i p . Further, it is noted that cot \pe lA 

The reflection is represented on a 
Poincar^" sphere in Figure A-5. The intersec- 
tion of the Si axis with the sphere (the 
point Ip , with Cartesian coordinates 1 , 0, 

0) represents the orientation of the x axis, or 
the plane of incidence, as I p corresponds 

to zero ellipticity x and zero azimuth <p by 
equations (A-31). The point I s has Carte- 
sian coordinates (-1, 0, 0) corresponding to 
zero ellipticity and an azimuth of -nj 2 ; 
therefore, it represents the orientation of 
the y axis, or the plane of the surface. As in 
the case of the double refracting plate, the 
phase change rotates the point P , which 
represents the polarization of the incident 
light, around the Si axis through an 
angle A , the phase change produced by 
reflection, to the point L , which has the 
coordinates 


r $ /r p and that tana.= a 2 /ai . 



Figure A-5. Representation of metallic reflec- 
tion on the Poincare sphere, side A [48] . 
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Si 1 


1 - tan 2 a 
1 + tan 2 a 


S 2 ' 


2 tan a cos (5 + A) 
1 + tan 2 a 


and 



(A-39) 


The point L is further translated to the point L' , whose coordinates are given in 
equations (A-38), by the relative absorption of the components caused by reflection. 


From equations (A-38) and (A-39), 
we have S 3 "/S 2 " = S 3 ' /S 2 ' = tan (5 + 
A) ;i.e., points L and L' lie on the great 
circle given by the locus of points lying on 
the Poincare sphere satisfying S 3 /S 2 = tan 
(5 + A) . This locus is shown in Figure A-6 
and is given by the intersection of the sphere 
with a plane normal to the S 3 axis and 
passing through the point L . Thus the 
point L is moved along the above great 
circle to the point L’ to correctly represent 
the polarization of the reflected light. Let us 
now derive the location of L' on the great 
circle. 

The coordinates of Ip and I $ are 

(1, 0, 0) and (-1, 0, 0), respectively, and 
those of L are given by equations (A-39). 
Therefore, the lengths of the 



chords IpL and I $ L are 


Figure A-6. Representation of metallic 
reflection on the Poincare sphere, side B [48] . 


)/ 1 - tan 2 a \ 2 [2 tan a cos (6 + A) 1 2 

jl 1 1 + tan 2 aj + [ 1 + tan 2 a J 

[" 2 tan a sin (6 + A) 1 2 ( 
+ \ 1 + tan 2 a j | 




4 ta n 4 a + 4 tan 2 a 
(l + tan 2 a ) 2 J 


2 tan a 


1 1 + tan 2 a) 


and 


I L 
s 


_ 1/ 

4 + 4 tan 2 a /2 _ 2 

(l + tan 2 a) 2 , 2 •. ‘/ 2 

L' ' J ( 1 + tan 11 a) 


(A-40) 


M220 
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so that • : 

. I L • , , _ . : . 

= tana . • ■ . (A-41) 

I i-t 

s 

IpLl s - is a right angle since it is inscribed in a semicircle; therefore, IpI s L = a . Note 
that tan a = a 2 /a t by equation (A-2). 

Now, the position of L' is determined by the angle IpI s L' • The tangent of this 

angle again gives the ratio of the amplitudes of the components of the reflected light, in 
analogy with equation (A-2). This ratio is given by tan a/tan \p . Therefore, 

I L» 

7^p- f = tan (i I L’) = tan a/tan i p . (A-42) 

1 .J-j p s 

s 

Using equations (A-38) as in equations (A-40) it is found that 


I L’ = 
P 


2 tan a cot ip 


1 / 

( 1 + tan 2 a cot 2 $) 2 


Moreover, using equations (A^tl) and (A-42), 
(I L )/(I s L) 

tan * - (yTyrnn ; 


(A43) 


(A44) 


in general, the ratio of the chords between a point representing the state of polarization of 
a given light wave and two diametrically opposed points on the sphere is equal to the ratio 
of the amplitudes of the components when the electric vector of the light is resolved along 
axes represented by these two opposed points. 

Let us now summarize the Poincare’ sphere representation of the reflection of light 

from a surface with known complex coefficient tan \p e 12 ^ . The point P , representing the 
polarization state of the incident light, is moved by the reflection process through an 
angle A on the sphere in a plane normal to the Si axis to the point L . The 
chords IpL and I S L are measured -and tan a calculated by equation (A-41). The 
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chord IpL' is calculated from equation (A-43), or the angle IpI s L' is calculated from 

equation (A-42), thus determining the state of polarization of the reflected light. It will 
later be shown how this representation can be used for the calculation of the effects 
occurring on reflection. First, however, still following the treatment of McCrackin, et al. 
[48] , the instrumentation will be described briefly. 


2. Instrumentation 

The various components of the ellipsometer are shown in Figure A-7. Collimated 
monochromatic light, that of the mercury green line (546.073 nm), is used. The polarizer, a 
Glan-Thompson prism mounted in a graduated circle, serves to polarize the source light. 
The compensator, also mounted on a graduated circle, is a mica quarter-wave plate; it is 
used to convert the linearly polarized light into elliptically polarized light. The incident 
light has an azimuthal angle and ellipticity predictable from the polarizer and compensator 
settings. In general, this azimuth and ellipticity will be changed by reflection from the 
sample. 


The aperture is an adjustable opening allowing a variation in the surface area 
examined and in the amount of light reaching the phototube. The analyzer, a second 
Glan-Thompson prism mounted on a graduated circle, can be rotated until a minimum 
intensity as indicated by the photometer is reached. If some ellipticity exists in the light 
reaching the analyzer, light extinction cannot be obtained by rotating the analyzer alone. 
In this case, the polarizer and analyzer must be adjusted alternately to remove this 
ellipticity and obtain extinction. Such an extinction, or null point, indicated that the light 
reflected from the sample is plane polarized. 

It should be noted that while alignment of the ellipsometer is not too critical for 
thickness and refractive index measurements on thin films, it is critical for measurements 
on surfaces. A step-by-step alignment procedure has been presented by McCrackin, et al. 
[48], 


3. Determination of Reflection Coefficient 

Tire reflection of light from a surface is characterized by the complex reflection 

coefficient tan i p e lA , and hence by the two quantities and A . Still following 
McCrackin, et al. [48] , it will later be shown how the refractive index and thickness of 
films on a surface are calculated from A and . Here, however, we shall deal with 
how A and are determined from readings on the analyzer (A) and polarizer (P) scales. 

For any given surface there is a multiplicity of polarizer, compensator, and analyzer / 
scale readings that produce extinction by the analyzer of the reflected light from the 
surface. In explaining these numerous readings and how A and i p may be computed 
from them, it should be remembered that all azimuthal angles are measured positive 
counterclockwise from the plane of incidence when looking into the light beam, and that 
the compensator is generally set so that its fast axis is at an azimuth of ±7r/4 . The present 
discussion is for an instrument such as that shown in Figure A-7, with the light passing 
through the compensator before reflecting off the sample. 
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LIGHT SOURCE & COLLIMATING LENSES 



1 


Figure A-7. The component parts of an ellipsometer [48]. 












The various readings fall into four sets or zones, two with the fast axis of the 
compensator set at -n/4 , numbered 1 and 3, and two with it set at -br/4 , numbered 2 and 
4. In each zone there is one independent set of P and A readings, or a total of four. 
However, since both polarizer and analyzer may be rotated by n without affecting the 
results, there appear 16 P and A settings, indicating a null point. Further, since the 
compensator may be rotated by n without affecting the results, there are finally 32 
possible sets of null point readings. 

In calculating A and i p , it is useful to calculate the three quantities p, ap , and 
a s from the P and A values, p being related to the P readings, ap being related to 
the A readings in zones 1 and 4, and a s being related to the A readings in zones 2 and 3. 

These quantities and their relations will now be examined, following the Poincare sphere 
treatment of Winterbottom [50] . 

The following definitions and rules for use of the Poincare sphere, which was shown 
in Figures A-3 through A-5, should be briefly noted: 

a. The Poincare sphere is a sphere of unit radius. 

b. Use of the sphere may be compared to use of a globe in navigation. The 
sphere has north and south poles (which lie on the S 3 axis) and is divided by latitudes, 
longitudes, great circles, and other arcs as needed. 

c. . Latitudes on the Poincare 
sphere show the ellipticity of the light or the 
ratio of the minor (b) axis to the major (a) 
axis of the ellipse. By definition, x is the 
ellipticity. Then , 2x is the latitude on the 
sphere and tan x = b/a . For example, circu- 
larly polarized light has b/a = 1 . Then tan x 
= 1 and x = 45 deg . The latitude for circu- 
larly polarized light is then 2x = 90 
deg , which is at the poles. For linearly 
polarized light, b/a = 0. Then, tan 
X and x are zero, and the points repre- 
senting this polarization lie on the equator. 

Elliptically polarized light has x between 
zero and one; therefore, all points except 
those on the equator and at the poles repre- 
sent elliptically polarized light (Fig. A-8). 

d. Longitudes represent the 
angle of inclination of the major axis of the 
ellipse (/> . At 0-deg longitude, which is the 
point Ip in Figure A-9 [51], the major axis 

is parallel to the equator. At 180 deg, or Figure A-8, Poincare sphere, the general 

point l s in Figure A-9, the major axis is significance of the different areas. 

perpendicular to the equator, <p = 90 deg . 

Longitudes are measured clockwise from 0 to 360 degj'thus, the angle of the major axis 
turns from 0 to 180 deg. 
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, e. The sense of rotation of the ellipse is counterclockwise in the southern 

hemisphere and clockwise in the northern hemisphere, using the traditional sense. 

The compensation process will now be given in terms of the Poincare 7 sphere 
representation [51]. 


In Figure A-9, as in the case of metallic reflection, Ip marks the plane incidence 
and I s marks the plane of the surface. In the ellipsometer, the corresponding planes are 

horizontal and vertical. Point Cf marks the 
fast axis of the compensator, which is set in 


Figure A-9. Poincare sphere, 
the ellipsometer compensation 



this case at 270/2 or 135 deg to the plane 
of incidence, which is the same as at an 
azimuth of -tt/ 4 . The fast. axis of the com- 
pensator is its optical axis. Point C s marks 

the slow axis of the compensator, which is 
180/2 or 90 deg from the fast axis. Point P 
is the polarizer angular reading, which in 
this case is roughly 40/2 or 20 deg clock- 
wise from Ip . 

The phase change caused by the 
compensator, 5 C , changes the linear polar- 
ization form at P into the elliptical form 
at Z . The phase change 5 C is 7r/2 for a 

perfect quarter-wave plate. As explained in 
the previous treatment of the double 
refracting plate, such a phase change rotates 
polarization forms clockwise on the sphere 
about the Cf - C s axis as one looks 

from Cf - C s through an angle equal in 


process [51]. magnitude to the phase change. 


Reflection of the elliptically polarized light Z from a metal surface produces a 
phase change A which rotates Z to L about the I n : L axis. In addition to the phase 

, . - v s ■ ; > .■ . 

change A , a relative absorption of the E x and Ey components occurs which shifts the 

linear polarization from L to L' . This L' form is extinguished by setting the 
transmission plane of the analyzer A at 180/2 or 90 deg from.- L' . In this -case, A. is 
roughly 60/2 or 30 deg. ' / >• • ‘ ! • •• v ; : 


As noted before, polarizing prisms and quarter- wave compensators can be rotated 
by 7r without affecting the polarization forms; To: avoid confusion, ellipsometer readings 
are grouped into four sets or zones. ‘ .'A. 

.Figure A-10 shows plane views of the Poincare; sphere which illustrate; how the 
four zones arise. Zone 1 corresponds to the diagram in Figure A-9. The P. and A readings 
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ZONE 1,C f ATT35 de§ ■' 


ZONE 2. C f AT 45 deg 





in Zone 1 will be designated p and ap , respectively. Since P and A can be rotated 

by t t and still give the same polarization form as before, there are four possible readings in 
this zone. These readings are combinations of P = p, p + n and A = ap , ap + ir .In Zone 

2, the case shown gives P = 7r/2 - p, 3rr /2 - p, , and A = a g , a s + v . Zone 3 has P = p + 
7r/2 , p + 37r/2 and A = n - a s , 2 jt - a s . Zone 4 has P = 7r-p,27r-p and A = rr - ap , 2n 

-v 

Zone T. The fast axis of the compensator is at -7r/4. The polarizer plane of 
transmission makes an angle of +p with the plane of incidence. The analyzer plane of 
transmission makes an angle of +ap with the plane of incidence. 

Zone 2. The fast axis of the compensator is at -Hr/4. The polarizer plane of 
transmission makes an angle of n/2-p with the plane of incidence. The analyzer plane of 
transmission makes an angle of +a g with the plane of incidence. 

Zone 3. The fast axis of the compensator is at -tt/ 4, The polarizer plane of 
transmission makes an angle of p + rr /2 with the plane of incidence. The analyzer plane of 
transmission makes an angle of -a s with the plane of incidence. 

Zone 4. The fast axis of the compensator is at +n/4. The polarizer plane of 
transmission makes an angle of -p with the plane of incidence. The analyzer plane of 
transmission makes an angle of -ap with the plane of incidence. 

Thus, the term p is the angle between the polarizer plane of transmission and 
either the plane of incidence or the plane of the surface, while a p and a s represent the 

angle between the analyzer plane of transmission and the plane of incidence, being labeled 
either a p or a s , according to whether p is the angle between the polarizer plane of 

transmission and the plane of incidence or the plane of the surface, respectively. The 
relations of P-and A to p and ap or a s are summarized in Table A-l . 

We shall now relate p , a„ , and a c to A and . 

P 5 .. 

Winterbottom [50] has solved the spherical triangles involved in the Poincare 
representation of the compensation process, as shown in Figure A-7, and has obtained 

,-tan A = sin 5 cot 2p 

and - (A-45) 

tan 2 ip = tah a tana , 

p S ( » 
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TABLE A- 1. RELATION OF P AND A READINGS 
TO p, a p , AND a s [48] 


Zone Compensator 


“7t/4 

P 

3 P 


P + 77 

' 

3 P 


P 

a p + 77 


P + 77 

a p + 77 

-77/4 

P + ff/2 

77 -a s 


p + 3^/2 

77 -a s 


p + rr/2 

277 -a s 


p + 37r/2 

■ 

277 -a s 

+7r/4 

7r/2 -p 

. a s 


3tt/2-p 

a s 


rr/2-p 

a S + 77 


3ir/2 - p ; 

3 s + 77 

■hr/4 

7T -p 

77 -a p 


27? - 77 

77 -3 p 


77 - p 

277 - a p 


277 - p 

2t 7 - a p 






where 5 is the relative retardation of the compensator. For a perfect quarter-wave 
plate (5 = rr/2), these reduce to 

A = 7r/ 2 + 2p 

and - (A-46) 


ib = a = a .... ...... 

p s 

Thus, by using Table A-l and equations (A-45) or (A-46), it is possible to 
compute p and ap or a $ from P and A , and hence compute A and <// . 

In practice, the values of p , a p and a g found from the P and A readings in the 

various zones are not identical. Experience shows that the averages from Zones 1 and 3 
check very closely with the averages from Zones 2 and 4 for both p and a , so 
measurements are usually taken in Zones 1 and 3 only and averaged. 


4. Computation of Film Indices and Thicknesses 

A typical thin Film under study consists of a film of index n 2 and 
thickness d 2 on substrate of index n 3 immersed in a medium of index nj , as shown in 
Figure A-l 1. Let all media be isotropic and n, represent a real index of refraction, while 



Figure A-ll. Reflection and refraction of light by one thin film. 
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n 2 and n 3 may be complex. This system is treated in Appendix B, where it is shown 
that the parallel and normal reflection coefficients for light incident at the boundary 
between the immersion medium and film are 



n 2 cos - nj cos 0 2 
n 2 cos <pi + cos 0 2 


and 


(A-47) 


s 

r 12 = 


nj cos (pi - n 2 cos 0 2 
n^ cos + n 2 cos 0 2 


Similar expressions hold for the reflection coefficients at the boundary between film and 
substrate, rP 3 and r| 3 . The cos 0j values heeded for these reflection coefficients are 
given by Snell’s law; i.e., 


cos 0 2 = [1 - (sin 0! n^n^ 2 ] 


% 


(A-48) 


etc. 

It is also shown in Appendix B that the total reflection coefficients that include the 
contributions of reflections from lower boundaries are given by: 

R P = + r 23 exp D 

P p 

1 + r 12 r 23 ex P D 

and (A-49) 

R s _ r 12 + rf 3 exp D 

g g * 

1 + r 12 r 23 ex P D 

where D is -4rr i n 2 cos 0 2 d 2 /X . 

The ratio of the total reflection coefficients defines p , where 

p = R P /R S (A-50) 
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which may be further expressed as 


P 


tan 



(A-5 1 ) 


As has been shown, we can calculate p from the ellipsometer readings. 

Therefore, by knowing n! and n 3 and by making ellipsometer measurements, we 
can determine refractive indices for films of known thicknesses or thicknesses of films of 
known indices. Note that the equations cannot be solved for d and n 2 in closed form. 
The normal procedure is to assume a series of refractive indices and to calculate thicknesses 
from the experimental measurements. Since manual calculations of this type are 
impractical because of the complexity of the equations, the problem is normally solved by 
computer techniques. In this study, the ellipsometer computer program of McCrackin [36] 
was adopted to Fortran 4H by Mr. J.L. Zurasky, MSFC, for use in the MSFC XDS Sigma 
Five Computer. This program calculated the refractive indices used in this study from 
ellipsometric measurements made on a Gaertner Model LI 1 9 ellipsometer. 
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APPENDIX B. DERIVATION OF THE DRUDE EQUATIONS 


Let us first derive Fresnel’s formulas for the reflection and refraction of light on a 
boundary between two isotropic dielectrics. We proceed as follows. 

Consider a plane light wave incident on a reflecting boundary. The medium from 
which the light ray is incident has permittivity e 0 and refractive index n 0 , and the other 
medium, in which the light ray is refracted, has permittivity e, apd refractive index n, . 
The permeability of both mediums is that of a vacuum, m v • Vector orientation is as shown 
in Figure B-l . 


\ 



Figure B-l . Reflection and refraction of light at a plane boundary 
. between two dielectrics. 

It should be noted that the velocity of light in the first medium (with index 0) is 
given by 

V 0 = 1/VerfT (B-D 
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and in the second medium (with index 1 ) by 


V t = l/ve^T 


(B-2) 


By using the complex form, where 


— — lx 

E = electric field strength = E 0 e 


(B-3) 


in which 






: 0,z 


= Ae^ 1 


Be ifi 2 , 


Ce i5 a , 


= oo(t - z/v), where z is the direction of propagation , 


5 , , 5 2 , 5 3 = phases , 


r time , 


= velocity of light in medium investigated , 


co 


GO 


angular frequency , 
2 irf , 


f = frequency , 

the vector plane wave of the incident light can be written as 


Eq = A 0 e 1T ° 


(B-4) 
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the vector reflected wave as 




(B-5) 


and the vector refracted wave as 


E? = A,e lT ' 


Note that for equations (fM) through (B-6) 


d / xsin0 o + zcos0 o \ 
t 0 - cult - -I , 

' Vrt ' 


r / xsin05 + zcos0j\ 

T 0 _ w it ~ ) » 

X V Q / 


and 


d 

T j — CU 




xsin^j + zcos^j 


) ■ 


(B-6) 


(B-7) 


(B-8) 


(B-9) 


where 4 > 0 , <Po \ and 4> i are defined as shown in Figure B-l . 

The amplitude vectors in equations (B-4) through (B-6) are divided into 
components normal to and in the plane of incidence; i.e., Ag. is divided 

into Aq p and Aq s , Rq is divided into Rg p and Rg s , and Aj is divided 
into A, p and A 1$ . 

Now, the above vector light waves will be written out into components onto the 
coordinate axis, noting that the components of the s amplitudes lie in the y direction and 
that the p components have projections onto the x and z axes. This may be seen from 
Figure B- 1 . 
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Thus, the following equations are obtained for incident light 

. d 

pd . 1T 0 

Vx = A o fP cos V 


F d _ ' 1T 0 

0,y 0,s e 


and 

. d 

d 1T 0 

E 0,z = - A o,p Bm V 

For reflected light we have 


r 1T 0 

E. = R cos 
0,x 0,p y 0 


< - E 0 
0,y 0, s 


and 


0, z 


-R. sin^'e 

0,p y 0 


. r 
lr o 


and for refracted light, 


h A lT i 

E l,x = A IlP COS V 

E d = A e" 1 , 

i»y i»s 


and 


d A . lT l 

E. - -A, suKpe 
l,z l,p r l 


CB-10) 
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Now, by Maxwell’s equations for isotropic homogeneous dielectrics containing no electric 
charges, which are 


and 




9H 

9E 

9E 

X 

= -Z - 

z 

v 9t 

9z 

3y 

9H 

9E 

9E 

y 

z 

X 

V 9t 

9x 

9z 

9H 

9E 

9E 

z 

X 

y 

v 9t” 

3y 

9x 


(B-l 1) 


and by integration, using equations (B-10) in equations (B-l 1), the following equations are 
found for incident light 


it. 


d , r— 0 

H = -n ve [i. A cos <pc 
0, x 0 v v 0, s r 0 


and 


H = n /\i A e 0 • 

0,y 0 v v 0,p 


. d 

d , r~ 1T n 

H n = n 've /M A. sin^ e u 
0,z 0 v v 0,s 0 


For reflected light we have 


1T„ 


H r = -n 'Ve /n B. n cosb’e 0 , 

0,x o v v 0,s 0 


. r 

1T„ 


= n n/7 /h e 0 , 

0,y 0 v v 0,p 


and 


. . r 

1T„ 


H* = n Ve /(x R sin^’e 0 , 

0,z 0 v v 0,s 0 ’ 


(B-l 2) 
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and for refracted light, 


. d 

IT, 


d IT 

H, = -n N/e /u A cos0 e x 
1 , X 1 V V l,s Y 1 


IT. 


H, = n, 's/e 7/ u A, e 1 
l,y 1 v v l.p 


and 


H 


d r—r a • 1T l 

, - n N€ /u A sm<6 e 

l,z 1 v"v 1, s *1 


Now, for two dielectrics, the tangential components of E and H and the 
components of B and D are continuous across the interface [52] . Therefore, 
write 


d r d 

E + E = E 
0,x 0,x l,x 


d : r d 

H a + H = H , 

0,x 0.x ,1.x 


d _r _d 
E + E = E 

O.y o.y l,y 


H d + H r = H d 
o.y o,y i,y 


e (V* + E* ) = e E d , n (n d + H r ) = p H d 

OV O.z 0,z/ 1 1, z v\ O.z O.z/ v l,z 

Now, substituting equations (B-l 1 ) and (B-l 2) into equations (B-l 3), we obtain 


. d . r . d 

1T 0 lr 0 1T 1 

A cos^ e + R costfde = A" cos <p e 
0,p 0 o.p 0 l.p 1 


. d . r . d 

lr o lr 0 lT i 

A e u + R e u = A e 1 , 

0,s 0, s 1, s 


j 

i T u iT. ir^ 

(A„ sind e ® + R„ sintf/e u ) = e A sinci e , 

V o.p r 0 0,p r 0 / 1 l,p r l 


. d 
it. 


d 

IT. 


IT IT IT 

n (A cos 0 e ^ + R cosd’e ^ = n A cos <b e , 
0\ 0, s 0 0,s 0/1 l,s 1 


normal 
we may 

(B-l 3) 

(B-l 4) 
(B-l 5) 
(B-l 6) 
(B-l 7) 
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(B-18) 


. d 

XT, 


. r 
it. 


. d 

IT . 


n (a e ® + R e = n A e * , 

0\ 0,p 0,p / 1 l,p 

. d . r 

IT 1T^ 

n (A.„ sin0 e • + R sintf’e V= n A, sind e 

(A 0, s r 0 0,s r 0 / 1 l,s r l 


. d 

IT. 


(B-19) 


Now, by the law of reflection, it is known that- cos 4> 0 ’ = -cos 4> 0 and sin <p 0 - sin 
<t> o' ,and by Snell’s law of refraction it is known that sin 0 o /sin 0, = n,/n 0 . Therefore, 
equation (B-19) may be written as 


IT, 


IT 


XT. 


A e 
0,s 


R e 
0,s . 


A e 

l.s 


which is the same as equation (B- 15). Noting that e 0 
(B-l 6) may be rewritten as .. ... 


= n 0 2 and ex = ni 2 , equation 


n n (a 
0\ 0 ,] 


. d ", ,r 

1T o p 1T o 

3 + R„ e 

. Q.P 




n A e 
1 l,p 


d 

1T 1 


which is the same as equation (B-18). 


• " Thus" four conditions remain that must be satisfied at z ;, = 0 for all t . These 

requirements are satisfied by the following conditions at , z - 0 . . s 

a. •' Phase condition : •• * ••• •• : 


T ' — T —T 
0 0 1 


Amplitude conditions '•••• ••• - 1 

( A - R„ ) cos0 = A, cos6 
v 0,p : . 0, p r () l,p / 1 

c V - _ ! / 9 ' 

(A + IT?-' ) n ’. = A h • , 

0,p 0,p y 0 l,p l • 


+ R = A 


0,s 0,s l,s 

( V'V )n o 0 ‘**o* A i,s n i cos ' t i 


(B-20) 


(B-21) 
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and for the reflected wave 


E*" = R exp 
0 0 H 


ioj 


(*- 


x sin0^ - 


z cos 0 


) 


IT 
^ . 0 


- V 


(B-27) 


In the film there are again two waves, the wave refracted by the first boundary, 
which includes all waves passing in the same direction after reflection at z = 0 and z = d , 
and the waves reflected at z = d . The wave refracted at z = 0 is written as 


-d - 

E = A exp 


x sin0„ + z cos 0, 

“ ( l r— — ) 


. d 

a lTl 
“ A l e 


(B-28) 


and the wave reflected by z = d is written as 


= R^exp 


xsin0 - z cos 0 

“(t- — ) 


IT, 




(B-29) 


In the third medium, we have all waves refracted by z = d . This set of waves is 
regarded as one wave which is written as 


"d 7 
E 2 = A 2 6XP 


ioj (t - 


xsin0 + z cos0 


IT, 


= V 


(B-30) 


In general, Aq , Rq , Aj , Rj , A 2 , and R 2 are complex numbers. 

Now, at boundaries I and II, we may apply Fresnel’s formulas, equations (B-22) 
through (B-25), which yield the following relations. 


I 1 tan (0 - 0 ) 

r » = R /A 1 = 

p 0 ,p 0 , p • tan (0 + 0 ) 


r’ = R /A - = - — — T~ 

s 0, s 0, s sin(0 c + 0 ) 




(B-31) 
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t» = A 1 / A 1 = - - 

p l,p 0,p sin(0 o +0 i )cos(0 o - 0 ) 

2sln* cos* 

t = A /A = 7 r 

s l,s 0,s sin(V> 0 + 0^ 

II II 

r M /A U = — . 1 \ 

P l,p l,p tan(0 i + 0 2 ) 

r" - R/ /A = - 7-^ 

s l,s l,s sin(0 + <p ) 

J- u 


^ sin © cos o 
1 r 0 


(B-32) 


(B-33) 


t" - /A^ = — i 

p 2,p l,p sin(0 + 0 2 ) cos (0^ - 0 2 ) 

2sin0 cos 0 

t" = A„ /A 2 1 


2sin0 cos0 1 


(B-34) 


s " J 2,s'"l,s sin(0 i + 0 2 ) 

Note that equations (B-10) are still valid for the components 

of E 0 d , E 0 r , and Ej d . Analogously, we may write for the components 

of E j r and E 2 T the following: 


. r 

IT 

r 1 

E , = -R, cos0 e 

1, x l,p r l 

. r 

r 1T 1 

E = R e 

l.y l,s 


E = -R sin0e 
l.z l,p 


. r 
1T 1 


(B-35) 


d 1T 2 

\x = A 2,p COS V 

. d 

d 1T 2 

E 2,y =A 2,s e ’ 
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. d 

d . 1T 2 

E 2,z = -\ P Sm V 


(B-35) 

(Con’t) 


Equations (B-12) for the components of Hq^ , Hq 1 , and Hj a are also still valid. 


H d 


Analogously, we may write for the components of H j r and the following: 


H = n N le /u. R cos 0 e 
1, X 1 v 1, s r l 


. r 
1T 1 


H*" = n N/e '/n R e 

l,y 1 v v l,p 


. r 
1T 1 


H, = n N / e /p R„ sin0 e 
1, Z 1 v v l,s 1 


. r 
1T 1 


IT, 


H = -n N/e/p A cos0 e 
2, x 2 v v 2,s r 2 

. d 

IT 

d , T~ , 2 

H = n Ve /p A„ e , 

2,y 2 v v 2, p 


H„ = n Ve/p A sin0 e 
2, z 2 v v 2,s r 2 


. d 
1T 2 


(B-36) 


As boundary conditions we have the continuity of the tangential components of 
the electric and magnetic field strength over boundaries I and II. This means that 
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f 


and 


E" 

+ 

E 

= E 

+ 

E 

0, x 


0,x 

1.x 


l.x , 

d 



d 


r 

E„ 

+ 

E 

= E, 

+ 

E, 

o.y 


o.y 

i.y 


i.y 

v d 


r 

d 


r 

H 

+ 

H 

= H 

+ 

H 

0,x 


0, x 

l.x 


l.x , 

v d 


„ r 

d 


r 

H 

+ 

H 

= H 

+ 

H 

o.y 


o.y 

i.y 


i.y 


E d + E r = E d 
l,x l,x 2,x 


d r d 

E + E* = E n 

I.y l.y 2,y 


d r d 

H + H = H 
l,x l,x 2,x 


H* + H* = H d 

i.y i.y 2,y 


(B-37) 


(B-38) 


Note that because of the discussion following equation (B-19), it is necessary to consider 
the requirements for continuity of the normal ' components of the electric and magnetic 
induction across the boundaries, since these conditions are contained implicitly in 
equations (B-37) and (B-38). 

The above boundary conditions are satisfied by the following assumptions: 

a. The phase conditions tq^ = tq = r and r = r j r = r are satisfied 

at z = 0 and z = d , respectively. 
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b. The following amplitude conditions are satisfied: 

At Boundary I (z - 0) 

(A 1 - R? ) cos </> = (A* - R 1 ).cos0 1 

V 0 , p 0,p ; r 0 v l,p l,p' y i 

(A 1 + R 1 ) n = (A 1 + R 1 )n 

v 0,p 0,p 0 l,p l,p 1 


A* + R n = A ^ + R \ 
0,s 0,s l,s 1, s 


(A. 1 - R 1 ) ri cos <p n = (A* - R„ ) n cosd 

v 0,s 0,s ; 0 r 0 v l,s l,s' 1 r l 

At Boundary II (z = d) 


II 

" R 1 r) 
l.P 

. A 11 

cos <6 = A n c 

l,p 

1 2,p 

II 

+ k i,p ) 

n = A n^ 

l,p 

1 2,p 2 

II 

+ H? 

= A I[ 

l,s 

l,s 

2,s 

II 

-< s > 

II 

n cos0, = A 

l,s 

1 r l 2, 


(B-39) 


(B-40) 


The amplitude conditions are obtained by substituting equations (B-10), (B-35), (B-12), 
and (B-36) into equations (B-37) and (B-38) and letting the phase conditions be fulfilled. 

Now, we must consider the phase relations between the amplitudes; i.e., we must 
introduce the principle of the interference of light into the calculation. From equations 
(B-26) through (B-30), it may be seen that we are concerned with a factor of co times the 
optical path difference, which is the product of the distance and the refractive index n(n = 
c/v, where c = speed of light in a vacuum)q.e., for a beam of monochromatic light that has 
been split and recombined, the phase difference between the two recombining beams 
is 2it/\ times (optical path difference), or 


Phase difference = 2n/X 




(B-41) 
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Now, Jenkins and White [54] show that the optical path in Figure B-2 is such that we have 
a phase difference of 


x = phase difference = ( 2 -tt/A. ) 2n^d cos 0^ 
Remembering that coz/v = (27r/X)cz/v = (2nj\ )nz , we have 


(B-42) 


l i.p 


= A 11 exp ( 
l,p \ 

,ku dcos<^ 

-A? e l/ * ix 
IfP 

) 

___ II / 

= R l, P e ^'( 

ko dcosri , 

_ 

= h“ e'‘ /jlx 
l.P 

V 1 ' 

- A 11 


-V 2 ix 

- A l,s e ' 

C 

s 


(B-43) 


Now we may combine equations (B-39) and (B-40) as: 

(A* -b' )co^„ =(a" 

v 0, p 0,p y r 0 ' l,p l,p ' v l 


/a 1 d 1 \ /.II Voix „II _1 /?ix\ 

(A + R \n = ( A e 1 + R e 1 )n 

V 0,p Oip/ 0 ' l,p -■ l,p '*1 


F J* iK a. t? 11 o-Hix 


A* + R* = A “ e' 2 *“ + R“ e 
0,s 0,s 1, s l,s 


(B44) 


/.I \ , /.II Kix II -y 2 ix\ 

( A 0,s ’ R 0,s) n 0 OOSi4 0 = ( A l. s e - E l,s e ) n i°° S V 

II II II 

Now noting that from equations (B-33) that Rj p = r£Aj p and Rj s 
and substituting e' 1(/2 * x = z and e 1/2 * x - 1/z, we obtain 


II 

r s^l,s 
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( A 0,p - B 0,p) cos ^0 - 


1 - r"z TT 

n ll 

— A„ cos0„ 

z l,p 1 


( A + R ) n 
V 0, p 0,p/ 0 


1 + r " Z TT 

P A II. 

A n 

z l,p 1 


A 1 + R l 
0,s 0,s 


1 + r”z TT 

s .II 

z l,s ’ 


(A n - R* \ncos6 = 
V 0,s 0,s/ 0 0 


1 " r ” z TT 

s .II 

A n cosri 

z 1, s 1 1 


Considering first the last two equations of (B-45), we may write 

I I 2 

A„ + R (l + r M z )n cos 6 

0, s 0, s v s ' 0 r 0 


-a' 

0,s 0,s 


A" - R„ (l - r"z ) n eos0 

S 1 1 


Now, note that by equation (B-3 1 ), we have 


r’ = - 


sin(0 o - 0j) _ n o cos0 o - n^os^! 
sin (0 O + n o cos0 o + n 1 cos0 1 


Therefore, 


1 + r f 


js _ cos 0 O 


1 - r T n 1 cos0 1 
s 


and equation (B-46) becomes 


A n + R n 
0,s 0, s 

A I „I 

A . - R„ 

0,s 0,s 


(1 + r’J z )(l + rf) 
( 1 - r’J z 2 ) ( 1 - rp 


(B-45) 


(B46) 


(B-47) 


(B-48) 
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Rewriting, we have 


A* [ ( 1 — r"z 2 )(l - r’ ) - (l + r"z 2 )( 1 + r T )] 
0,s ' s s s s 


= R 1 [ - ( 1 + r" z 2 ) ( 1 + r f ) - ( 1 - r” z 2 ) ( 1 - r* 

0, s v s s' s " ' s 


or 


R 


0,s 


0,s 


i6 


r e 
s 


[(1 + r"z 2 )(l + r’) - (l - r''z 2 )(l - r* )] 

bo b b 

[(1 + r"z 2 )(l + r* ) + (l - r’’z 2 )(l - r* ) ] 

S b b b 


r M z + r* 


r* + r"e 


-IX 


1 + r M r* z 
s s 


1 + r* r M e 
s s 


-IX 


Now considering the first two equations of equations (B-45) we have 


II 2 

A + R ( 1 + r M z ) n, cos 

0,p 0,p p 1 0 

II 2 

A„ - R (l - r"z )n cos0 

0,p 0,p v p 0 1 


Note that by equation (B-3 1 ) we have 

, _ tan (0Q - 0 t ) _ n 1 cos0Q - npcos^! 
p tan (0 O + 0j) ni cos 0 O + ^008 0! 

Therefore, 

1 + r’ 

p _ n t cos 0Q 

1 - n 0 cos 0 j 


) ] (B-49) 


(B~50j 


(B-51) 


(B-52) 


(B-53) 
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and equation (B-5 1 ) becomes 


A 0,P + R 0,p 
A o,P - <P 


( 1 + r”z )( 1 + r' ) 
P ' P 

( 1 - r” z 2 ) ( 1 - r’ ) 
P P 


(B-54) 


Rewriting we have 


A* [ ( 1 - r” z 2 ) ( 1 - r’ ) - ( 1 + r M z 2 ) ( 1 + r' ) ] 

0,p p P P P 

= R* p [ - (1 + r p 'z 2 )( 1 + r p ) - (1 - r”z 2 )(l - r^) ] 


(B-55) 


or 


R 


0,p _ 


o,p 


r e 
P 


i6^ ( 1 + r"z 2 )( 1 + r r ) - ( 1 - r”z 2 ) ( 1 - r’ ) 

P P P P 

(1 + r” z 2 ) ( 1 + ) + ( 1 - r” z 2 ) ( 1 - ) 


(B-5 6) 


2r M z + 2r T 
P P_ 

2 + 2r”r’ z 2 
P P 


r* + r M e 
P P 


-IX 


1 + r’ r”e 
P P 


-IX 


We have thus obtained the exact Drude equations for the reflection coefficients for 
a film-covered surface in the form 


R = 


i6 

re P 
P 


r f + r”e 
_ P P 


-IX 


1 + r* r"e 
P P 


-IX 


i5 


-IX 


R = 


r e 
s 


1 + r T r” e 
s s 


-XX 


(B-5 7) 
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where x is (27f/X)(2n! d cos 0 j) and r § ' , r s " , Tp' , and and Tp" are the well-known 
Fresnel amplitude coefficients. 

It should be noted that the derivation just presented treated the case of a dielectric 
(nonabsorbing) media. We shall follow the common approach of extending our results to 
include the case of conducting (absorbing) media by replacing the real index of refraction 
by a complex quantity (n + ik, where k = 0 for a dielectric media). The extension proves to 
be approximately correct and is in accordance with a careful development of the Fresnel 
coefficients by Stratton [55] . This extension is also discussed by Nestall [56] . 

We note further that the classical Fresnel equations, while generally inadequate in 
metal optics, are applicable here since the frequency involved is below the plasma 
frequency, Wp . At cop and above, the effects of polarization waves must be considered. 

For most metals, u > p lies in the ultraviolet. This facet of the Fresnel equations is discussed 

by Melnyk and Harrison [57], who point out that when polarization fields become 
important the classical Fresnel equations no longer apply. However, for frequencies up to 
and including visible light frequencies, there are no measurable bulk charge densities in 
metals; hence no polarization fields. 
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